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SUMMARY 


This  report  describes  the  effort  completed  to: 

• Determine  the  parameters  that  most  significantly  influence  airplane  stopping  distance, 

• Formulate  a model  that  could  be  used  to  predict  airplane  braking  distance. 

• c’  -iuate  the  Diagonally  Braked  Vehicle  (DBV)  and  Mu-Meter  with  respect  to  the 
criteria  resulting  from  the  model, 

• Recommend  the  most  logical  approach  to  providing  this  prediction  capability  for 
USAF  and  FAA  usage. 

During  Task  I,  the  Boeing  Brake  Control  Simulator  war  used  to  identify  the  characteristics 
affecting  airplane  stopping  performance.  Five  aircraft  models  (727-200,  7S7  Advanced, 
747-200,  C-141  A,  and  F-4E)  were  the  subject  of  a sensitivity  study.  Major  airplane,  landing 
gear,  ai  d brake-wheel- tire  parameters  were  varied  systematically  from  arbitrary  baseline 
conditions.  This  technique  was  used  to  determine  the  parameters  that  influence  braking 
distance  by  at  least  2%. 

The  initial  list  of  parameters  was  reviewed  and  reduced  by  eliminating  those  involving  pilot 
technique  and  by  grouping  interrelated  terms.  The  resultant  list  is. 

• Peak  available  ground  friction  (q) 

• Drag  device  effectiveness  (Cl/Cq) 

• Brake  application  speed  (v) 

• Air  densi./  (p) 

• Engine  idle  thrust  (Fe> 

During  Task  II,  a prediction  model  has  been  developed  that  correlates  with  measured  test 
data  to  within  ±59!-.  The  model  consists  of  a prediction  equation  expressing  the  relationships 
between  four  dimensionless  groups  of  factors  (pi  terms)  needed  to  define  the  braking 
phenomenon.  The  equation  is  of  the  form: 

- - ^ ,C|  Co  . Cq 

" i '-  a v "3'  v "4; 

where  w j = sg/v^,  rro  = p,  S3  = C^ICq,  ir^  = pv6/ Fcg2,  and  Ca,  Cj,  Co,  and  C3  are 
determined  from  particular  aircraft  test  results. 

The  DBV  and  Mu-Meter  were  evaluated  using  the  developed  specification  criteria.  Both 
vehicles  perform  their  design  functions,  but  both  fail  to  meet  the  aircraft  criteria  1 , eloped 
during  this  study. 
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In  the  recommended  approach,  the  prediction  methodology  consists  of  a stopping  perform- 
ance equation;  its  applicability  depends  on  a meaningful  runway-, ‘ire  friction  coefficient  and 
correct  weather  information. 


vi 


TABLE  OF  CONTENTS 


Section  Title  Page 

Glossary ix 

List  of  Abbreviations  and  Symbols  x 

I.  Introduction I 

II.  Identification  of  Parameters  7 

1.  Airplane  Flight  Characteristics 7 

2.  Runway  and  Environmental  System 7 

3.  Landing  Gear  System  7 

III.  Stopping  Distance  Derivation  II 

1.  Air  and  Transition  Distance  II 

2.  Braking  Distance 13 

IV.  Sensitivity  Study  Test  Procedure  15 

V.  Sensitivity  Study  Test  Conditions 19 

1.  Stability  Studies 26 

2.  Performance  Studies 26 

3.  Hydraulic  System  Studies 26 

VI.  Test  Results  27 

1.  Baseline 27 

2.  Sensitivity  27 

a.  Airplane  Flight  Characteristics 26 

1.  Landing  Weight  29 

2.  Center  of  Gravity  29 

(a)  High  and  Low  Center  of  Gravity 29 

(b)  Forward  and  Aft  Center  of  Gravity 29 

3.  Brake  Application  Speed  29 

4.  Spoiler  or  Drag  Chute  Deployment  45 

5.  Spoiler  or  Drag  Chute  Effectiveness 45 

6.  Engine  Idle  Thrust 45 

7.  Brake  Pressure  Application  Rate 45 

8.  Delay  of  Brake  Pressure  Application  45 

9.  Metered  Pressure 45 

b.  Runway  and  Environmental  System 46 

1.  Wind 46 

2.  Hot  and  Cold  Days  46 

3.  Rough  Surface 46 


vU 


TABLE  OF  CONTENTS  (Continued) 


Section  Title  Page 

c.  Landing  Gear  Systems 46 

1.  Brake  Fade 46 

2.  Torque  Peaking  47 

3.  Torque  Response  Frequency 47 

4.  Torque  Gain  47 

5.  Mu-Slip  Curve  Variation  47 

6.  Worn  Tire 47 

7.  Strut  Frequency 47 

8.  Vertical  Stiffness  and  Vertical  Damping  48 

d.  Hydraulic  Systems 48 

1.  Line  Diameter  48 

2.  Line  Length  48 

3.  Restricted  Return  48 

4.  Brake  Volume  Changes 48 

e.  Composite  Variation  of  Parameters  49 

VII.  Parameter  Evaluation  Criteria  51 

1.  Perfc  nance  Indices 51 

a.  Airplane  Braking  Distance  51 

b.  Braking  Distance  Efficiency  51 

c.  Developed  Mu  Efficiency 52 

d.  Skid  Index  52 

e.  Cornering  Index  53 

f.  System  Stability 53 

2.  Parameter  Rating  System 53 

VIII.  Simulator-to-Airplane  Correlation 55 

1.  727  57 

2.  737  57 

3.  747  00 

4.  C- 1 4 1 A 60 

5.  F-4 66 

6.  Operational  Anomalies  66 

IX.  Determination  of  Significant  Parameters  73 

1.  Parameter  Ratings 73 

2.  Significant  Parameters 75 


TABLE  OF  CONTENTS  (Continued) 


Section  Title  Page 

“ X.  Selection  of  Pertinent  Parameters 77 

1.  Brake  Application  Velocity 77 

2.  Aerodynamic  Lift  and  Drag  77 

5 3.  Wind  Velocity  77 

4.  Center-of-Gravity  78 

5.  Other  Parameters  78 

XI.  Development  of  Prediction  Model 83 

1.  Component  Equations  83 

2.  Generalized  Functions  85 

XII.  Model-to-Simulator  Correlation  93 

XIII.  Model-to-Simulator  Qualitative  Comparison 97 

1 . Predicted  Mu  from  Flight  Test  Data 97 

2.  Wet-Runway  Analysis 97 

3.  Data  Validity 97 

XIV.  Determination  of  Vehicle  Criteria 101 

XV.  Evaluation  of  DBV  and  Mu-Meter 105 

XVI.  Conclusions  and  Recommendations 113 

1 . Conclusions 113 

2.  Recommendations 113 

References  1 1 fi 

Appendix  A-Calculation  of  Approach,  Flare,  and  Transition  Distances  119 

1 . Approach  Distance  119 

2.  Flare  Distance  119 

3.  Transition  Distance  1 20 

' Appendix  B -Boeing  Brake  Control  Simulation  123 

1.  Hardware 123 

2.  Airplane  Dynamics  124 

t 3.  Strut  Dynamics  124 

4.  Wheel  Dynamics 124 

5.  Brake  Torque  Simulation 124 

6.  Ground  Force  Simulation 124 

7.  Nonlinear  Relationships 124 


ix 


TABLE  OF  CONTENTS  (Concluded) 


Section  Title  Page 

a.  Wet  Runway  Curve  125 

b.  Touchdown  Profile  125 

c.  Nonlinear  Torque  Gain  Curves 125 

d.  Mu-Siip  Curve  Variations 125 

e.  Power  Spectral  Density  Curve  of  Rough  Runway  126 

8.  Recorded  Data 126 

9.  Test  Equipment  126 

Appendix  C-Dimensional  Analysis  Technique 137 

1.  Form  of  Dimensional  Equations 137 

2.  Development  of  Prediction  Equations 137 

3.  The  Buckingham  Pi  Theorem 140 

4.  Statement  of  Dimensional  Homogeneity  140 

5.  Determination  of  Pi  Terms 140 

Appendix  D-Formulation  of  Component  Equations 145 

1.  Polynomial  Regression  145 

2.  Multiple  Linear  Regression 146 

Appendix  E-Formulation  of  Generalized  Functions 149 

1.  Determination  of  Functions  149 

2.  Conditions  for  Function  to  be  a Product 149 

Appendix  F -Model  Theory  Analysis  for  Ground  Vehicles 155 

1.  Model  Theory  155 

2.  Model  Distortion 157 

a.  The  Prediction  Factor 158 

b.  Multiple  Distortion  159 

c.  Compensated  Distortion  159 


x 


gjjS^gggSE  " 


LIST  OF  ILLUSTRATIONS 

Fi*?re  Title  Page 

No. 

1 Basic  Program  Plan  2 

2 Basic  Program  Technical  Approach  3 

3 Landing  Segments 4 

4 Approach  Distance  vs.  Glide  Slope  12 

5 Flare  Distance  vs.  Approach  Speed  12 

6 Transition  Distance  vs.  Approach  Speed  vs.  Time  to  Brake  Application  12 

7 Transition  Distance  vs.  Approach  Speed  vs.  Aircraft  Deceleration 12 

8 Simulator  Block  Diagram 14 

9 Baseline  Braking  Distance  vs.  Peak  Available  Mu 28 

10  Braking  Distance  Test  Results  . . . . 30 

11  727  Dry  Runway  58 

12  737  Dry  Runway  59 

13  747  Dry  Runway  61 

14  747  Wet  Runway  62 

15  C-l 41  Wheel  Speed  Transducer  Simulation 63 

16  C-l 41  Dry  Runway  64 

17  C-l 41  Wet  Runway  65 

18  F-4  Dry  Runway 67 

19  F-4  Wet  Runway 68 

20  727  Flooded  Runway 70 

21  737  Wet  Runway  Early  Brake  Application 71 

22  Block  Diagram  for  Talk  11  Analysis 79 

23  Effect  of  Wind  at  7T2  = 0.6  80 

24  Effect  of  Wind  at  it  2 = 0.2  81 

25  Mu-Efficiency  Curves 95 

26  Mu-Velocity  Curve  for  Wet  Runway 99 

27  Recommended  Ground  Vehicle  Criteria  103 

28  Comparison  of  Tire  Sections 106 

29  Mu-Meter  Evaluation  Ill 

30  Recommended  Prediction  Methodology  114 

A-l  Air  and  Transition  Distances 121 

B-l  Block  Diagram  of  Boeing  Brake  Control  Simulator  129 

B-2  Wet  Runway  Curve  131 

B-3  Touchdown  Profile  131 

B-4  Nonlinear  Torque  Gain  Curves 132 

B-5  Percentage  of  Peak  Available  Mu  vs.  Percentage  of  Slip 133 

B-6  Reduced  Frequency  Power  Spectral  Density  Curve  of  Rough  Runway  Simulation  . 134 

D-l  Plots  of  wj,  vs.  7T2,  ^3,  and  7T4  147 

D-2  Equal.on  Flow  Chart  148 

F-l  Mu-Meter-to-Aircraft  Aerodynamic  Similarity 161 


xi 


LIST  OF  TABLES 


1 Parameters  Influencing  Stopping  Distance  8 

2 Definition  of  Simulator  Parameters If, 

3 Baseline  Parameter  Values  Used  in  Airplane  Simulation  18 

4 Test  Conditions  20 

5 Parameter  Change  Values 23 

6 Correlation  Tests 5f> 

7 Parameter  Rating  74 

8 Reduced  Significant  Parameters 75 

9 Significant  Parameters 77 

10  Pertinent  Independent  Variables 84 

1 1 Summary  of  Component  Kquations 86 

12  Test  of  Validity  for  the  Function  to  ne  a Product 90 

13  Test  of  Validity  for  Constant  Term 91 

14  Summary  of  Prediction  Kquations 92 

15  Summary  of  Percentage  Errors 94 

16  Qualitative  Comparison  with  Flight  Test  Data 98 

17  Comparison  of  Aircraft  and  Automotive  Tires 107 

18  DBV  Correlation  Requirements  109 

19  Mu-Meter  Correlation  Requirements  110 

B-l  Pressure  Transducer  Information  125 

C-l  Parameters  for  Braking  Stop  Distance  Model 144 

E-l  Calculations  of  Validity  for  the  Function  to  oe  a Product  153 

E-2  Calculations  of  Validity  for  Constant  Term  1 54 

F-l  Prediction  Factor  6 vs.  Distortion  factor  a for  DBV  Evaluation 162 

F-2  Mu-Meter  Data K,3 

F-3  Mu-Meter  Data  Analysis IP4 


GLOSSARY 


sir  distance  summation  of  approach  and  flare  distances 

approach  distance  distance  an  aircraft  travels  from  threshold  to  initiation  of  flare 

baseline  airplane  aircraft  in  a three-point  taxi  a'titude  (during  braking)  and  of  typical 
landing  weight,  approach  speed,  CG  location,  landing  flap  setting,  and  'ngine  thiest 

braking  distance  distance  an  aircraft  travels  from  brake  application  to  a lower  safe-tumoff 
speed 

contamination  foreign  substance  on  a runway 

flare  distance  distance  an  aircraft  travels  during  flare  segment 

flare  load  factor  normal  aircraft  acceleration  experienced  during  flare,  expressed  in  g’s 

glide  slope  angle  in  degrees  defining  aircraft  flight  path  during  approach 

height  above  threshold  aircraft  alt  tude  at  start  of  approach  segment 

parameter  rating  index  rank-order  relationship  of  parameters 

peak  available  mu  computer  input  defining  maximum  value  of  friction  available  between 
the  tire  and  ground  during  test  condition 

pi  term  dimensionless  term  for  modeling  theory 

slip  relationship  between  the  braked  wheel  speed  and  the  synchronous  wheel  (airplane) 
speed 

stopping  distance  summation  of  airplane  approach,  flare,  transition,  and  braking  distances 

transition  distance  distance  airplane  travels  from  time  of  initial  main  gear  touchdown  to 
brake  application 


i 

\i 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 


a 

aircraft  deceleration 

AC 

alternating  current 

CD 

coefficient  of  drag 

CG 

center  of  gravity 

CL 

coefficient  of  lift 

DC 

direct  current 

Fe 

engine  idle  thrust 

^eo 

engine  idle  thrust  at  zero  velocity 

g 

acceleration  caused  by  gravity 

H 

height  above  threshold 

KE 

change  in  idle  thrust  with  velocity 

N,n 

flare  load  factor 

PB 

brake  pressure 

PRI 

parameter  rating  index 

PSD 

power  spectral  density 

s,  SB,  XA 

braking  distance 

SA 

approach  distance 

SF 

flare  distance 

ST 

transition  distance 

tb 

brake  torque 

TA 

time  to  brake  application 

v,  VI 

brake  application  speed 

vw 

headwind  or  tailwind  velocity 

xiv 


t 


W airplane  landing  weight 

7 glide  slope 

Mi  mu  peak  available  friction  coefficient 

t) s  braking  distance  efficiency 

>r  pi  term 

p air  density 

o slip,  percentage  of  slip 

u)  wheel  velocity 


xv 


.'fyrryf  "wry — ^^■~T-r!7V!^rT> 


SECTION  I 

INTRODUCTION 

At  some  time  during  the  approach  to  an  airport,  an  airplane  pilot  must  decide  whether  or 
not  he  has  enough  runway  length  to  stop  safely.  To  make  this  decision,  he  must  consider  or 
anticipate  a variety  of  conditions.  Weather,  runway  conditions,  uirplane  touchdown 
dynamics,  and  braking  system  capability  have  a major  influence  on  the  stopping  distance  of 
the  aircraft.  Thus,  the  pilot  must  nave  information  not  only  on  his  own  aircraft's  landing 
characteristics  but  also  the  environment  in  which  he  is  lunding.  Because  a successful  landing 
is  vital  to  the  safety  and  economy  of  airplane  operation,  there  has  been  a continuing  effort 
by  the  FAA  und  USAF  to  provide  pilots  with  information  so  that  he  can  predict  the 
stopping  performance  of  his  aircraft  upon  landing. 

The  desire  to  provide  this  information  requires  a thorough  understanding  of  the  parameters 
that  determine  stopping  distance.  To  determine  what  technical  data  is  required  by  the  pilot, 
a joint  NASA,  FAA,  und  USAF  program  was  ini,  ated  and  the  Boeing  Commercial  Airplane 
Company  wus  awarded  USAF  Contract  F33657-74-C-01 29  to  identify,  evaluate,  and  rate 
the  parameters  affecting  airplane  stopping  performance.  The  busle  approach,  shown  in  Fig- 
ures 1 and  2,  to  meeting  this  objective  wus  to  list  ull  parameters  which  could  possibly 
Influence  stopping  performance.  Then,  with  the  aid  of  a digital  computer  and  an  analog- 
hardware  brake  contiol  simulation,  the  parameters  were  varied  about  their  baseline  to 
determine  their  uctuul  effect  on  stopping  distance.  To  druw  some  general  conclusions  re- 
garding the  significance  of  individual  parameters,  five  uirplunes  were  studies  (727-200,  737 
Advanced,  747-200,  C-I41A,  and  F-4U).  The  process  of  identifying,  evaluating,  und  rating 
the  fuctors  as  described  above  lias  been  termed  a sensitivity  study. 

To  produce  meaningful  resuits  in  a sensitivity  study,  all  influencing  fuctors  must  be  con- 
sidered. Thus,  the  unulysis  of  stopping  performance  should  include  parameters  that  affect 
the  uirplane  while  it  is  on  the  runwuy  and  also  factors  that  determine  its  initiul  condition  al 
brake  application.  Thus,  runwuy  performance  involves  the  totul  uirplane-runwuy-weuther 
system  from  the  point  at  which  the  airplane  crosses  the  runway  threshold  to  the  point 
where  it  actually  stops.  As  a result,  the  stopping  distance  has  been  divided  into  four  distinct 
segments:  approach,  flare,  transition,  and  braking  (Figure  3).  The  approach  and  flare 
segments  are  air  modes  that  determine  the  dynamic  conditions  of  the  airplane  at 
touchdown.  Transition  involves  the  distance  the  airplane  travels  as  it  changes  from  an  air 
mode  to  the  braking  mode.  Braking  is  the  distance  during  which  the  airplane's  forward 
kinetic  energy  is  dissiputed, 

Determination  of  Federal  Air  Regulations  (FAR)  field  length  requires  an  accounting  of  each 
segment.  During  the  air  and  transition  modes,  the  distance  traveled  is  determined  largely  by 
airplane  flight  characteristics,  pilot  technique,  and  the  rules  governing  operation  during 
these  segments.  The  same  factors  influence  the  braking  distance  because  they  contribute  to 
the  airplane  dynamics  at  brake  application, 


Figure  1. —Basic  Program  Ptan 


at  Traction 
I data 


Figure  2.  -Basic  Program  Technical  Approach 


Distances  from  the  threshold  to  touchdown  can  be  readily  calculated  by  considering  the 
individual  segments  of  approach  and  flare,  but  a meaningful  braking  distance  can  be  deter- 
mined only  with  an  airplane  test  or  a brake  control  simulation.  Because  airplane  testing  is 
expensive  and  time-consuming,  the  simulator  is  a useful  tool,  both  to  evaluate  different 
antiskid  systems  and  to  analyze  sensitivity  where  each  parameter  must  be  closely  controlled. 
Results  of  a series  of  simulator  tests  can  be  applied  to  a dimensional  analysis  because  only 
one  influencing  parameter  is  varied  at  a time.  The  goal  of  the  dimensional  analysis  is  the 
development  of  a braking  distance  prediction  equation.  This  equation  must  include  all 
variables  having  a significant  influence  on  bilking  distance  during  the  normal  operation  of 
military  and  commerciat  lrcraft. 

Statistical  curve  fitting  techniques  can  be  applied  to  develop  a relationship  (or  equation) 
between  the  dependent  and  any  other  independent  variable  while  all  other  parameters  are 
held  constant.  This  process  is  repealed  for  each  variable.  These  component  equations  can 
then  be  combined  to  form  a prediction  equation,  provided  the  necessary  and  sufficient 
conditions  generated  during  the  analysis  are  met. 

The  Diagonally  Braked  Vehicle  (DBV)  and  the  Mu-Meter  are  two  ground  vehicles  in  use 
today  to  evaluate  runway  friction  conditions.  Various  methods  have  been  explored  and 
applied  in  an  attempt  to  correlate  the  measured  ground  vehicle  data  to  airplane  braking 
distance.  With  the  use  of  a dimensional  prediction  equation,  the  vehicles  can  be  compared  to 
the  dimensional  groups.  Vehicle  usefulness  or  modification  to  satisfy  the  dimensional 
prediction  equation  can  be  assessed. 
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SECTION  II 

IDENTIFICATION  OF  PARAMETERS 


Accurate  prediction  of  aircraft  stopping  distance  requires  that  all  significant  parameters  be 
taken  into  account.  The  many  parameters  determining  aircraft  stopping  distance  have  been 
divided  into  groups  based  on  the  nature  of  the  parameter  or  the  system  the  parameter 
affects.  The  three  major  groups  identified  are: 

• Airplane  flight  characteristics 

• Runway  and  environmental  system 

o Landing  gear  system 

A detailed  list  of  parameters  associated  with  these  groups  is  given  in  Table  1 . The  following 
brief  descriptions  of  the  major  groups  explain  the  significance  of  some  of  the  parameters 
listed  in  Table  I. 


1.  AIRPLANE  FLIGHT  CHARACTERISTICS 

The  airplane  system  is  the  mass  being  stopped.  As  such,  it  provides  the  Inertial,  propulsive, 
and  aerodynamic  forces  that  determine  the  amount  of  energy  to  be  dissipated  by  the  brakes. 
The  airplane  also  affects  braking  in  that  tire  loading  is  determined  from  airplane  geometry, 
dynamics,  and  aerodynamics.  The  pilot  al'fec'x  stopping  by  determining  the  initial  condi- 
tions of  the  aircraft. 


2.  RUNWAY  AND  ENVIRONMENTAL  SYSTEM 

The  runway  and  environmental  system  consists  of  effects  that  influence  the  airplane 
stopping  performance  but  are  external  to  the  aircraft,  These  include  atmospheric  conditions, 
altitude,  temperature,  runway  slope,  crown,  surface  roughness  (micro-  and  macro-texture), 
material,  and  contamination. 


3.  LANDING  GEAR  SYSTEM 

The  landing  gear  system,  as  it  affects  stopping,  consists  of  the  structure,  shock  strut,  and 
braking  system.  The  structure  transmits  the  forces  developed  between  the  tire  and  runway 
to  the  airplane,  while  the  shock  strut  is  the  primary  vertical  energy  absorption  system  on  the 
airplane.  The  flexibility  and  damping  characteristics  of  the  structure  and  the  dynamic 
characteristics  of  the  shock  strut  have  a pronounced  effect  on  stopping  system  performance. 
The  uraking  system,  consisting  of  the  wheels,  brakes,  tires,  control  system,  and  brake 
hydraulics,  is  of  primary  importance  in  determining  the  ground  distance.  Changes  to  the 
braking  system  will  vary  the  rate  and  magnitude  of  k inetic  energy  transfer 
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Table  1 .—Parameters  Influencing  Stopping  Distance 


I.  AIRPLANE  FLIGHT  CHARACTERISTICS 

1.  Aerodynamics 

i.  Brako  application  spaed 

b.  Aerodynamic  coefficient!  (lift  and  drag) 

c.  Spoiler  deployment  rate  and  effectiveness 

d.  Engine  idle  and  spindown  thrust 

2.  Geometry 

a.  Center  of -gravity  location 

b.  Landinq  gear  placement 

c.  Airplane  weight 

d.  Airplane  mass  moment  of  inertia 

e.  Wing  stiffness 

3.  Operational  characteristics 

a.  Route  structure 

b.  Pilot  braking  technique 

II.  RUNWAY  AND  ENVIRONMENTAL  SYSTEM 

4.  Runway 

a.  Roughness  (micro-  and  macro-texture) 

b.  Contamination 

c.  Slope 

d.  Crown 

e.  Mu 

5.  Atmospheric 

a.  Ambient  temperature 

b.  Pressure  altitude 

c.  Wind 

III.  LANDING  GEAR  SYSTEM 

A.  Brake  System 

- Brake  Actuation  System 
0.  Mechanisms 

a.  Pedal  crank  mechanical  advantage 

b.  Actuation  mechanism  mechanical  advantage 

c.  Mechanism  detailed  design 

d.  Cable  length 

e.  Cable  stiffness 

f.  Cable  preload 
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Table  1 .-Parameters  Influencing  Stopping  Distance  (Continued) 

III.  LANDING  GEAR  SYSTEM  (Continued) 

7.  Hydraulics 

a.  Hydraulic  pump  capacity 

b.  Pump  recovery  rate 

c.  Accumulator  pressure  volume  characteristic* 

d.  Fluid  bulk  modulus 

e.  Fluid  viscosity 

f.  Fluid  density 

g.  Fluid  temperature  coefficients 

h.  Brake  line  length 

i.  Brake  line  stiffness 

j.  Brake  line  diameter 

k.  Flow  restrictions 

l.  Hydraulic  fuse  type 

m.  Bleeding  characteristics 

n.  Return  line  steady-state  back  pressure- 

o.  Return  line  transient  back  pressure 

p.  Metering  valve  flow  capacity 

q.  Metering  valve  pressure  gain 

- Antiskid  System 

8.  Wheel  speed  sensor 

a.  Eccentricity 

b.  Drive  mechanism 

c.  Stator  mounting 

d.  Electrical  characteristics 

e.  Drive  angle 

f.  Demodulator  characteristics 

g.  Sensor  type 

h.  Sensor  detailed  design 


a.  System  type 

b.  Detailed  design 

10.  Servo  valve 

a.  Valve  type 

b.  First-stage  type 

c.  Flow  gain 

d.  Detailed  design 

e.  Lapping 

- Wheel,  Brake,  and  Tire  System 

11.  Wheel  and  Brake 

a.  Pressure  volume  characteristics 


a 
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Table  1 .—Parameters  Influencing  Stopping  Distance  f Concluded ) 


III.  LANDING  GEAR  SYSTEM  (Continued! 

b.  Lining  characteristics 

c.  Brake  detailed  design 

d.  Heat  sink  loading 

e.  Heat  transfer  characteristics 

f.  Retractor  spring  pressure 

g.  Wear  adjuster s 

h.  Brake  mounting  system 

i.  Wheel  design 

12.  Tire 

a.  Footprint  area 

b.  Inflation  pressure 

c.  Fore  and  aft  stiffness 

d.  Vertical  stiffness 

e.  Lateral  stiffness 

f.  Mass 

g.  Inertia 

h.  Tread  pattern 

I.  Thermal  properties 

j.  Rolling  resistance 

k.  Diameter 

l.  Damping  characteristics 

m.  Wear 

B.  Shock  Strut  and  Structure 

1 3.  Strut 

a.  Fore  and  aft  stiffness 

] b.  Strut  mass 

c.  Strut  effective  mass 

d.  Fore  and  aft  damping 

' e.  Vertical  stiffness 

f.  Metering  pin  design 

g.  Torsional  stiffness 

h.  Torsional  damping 

i.  Tolerance* 

j.  Strut  geometry 

14.  Truck 

a Mass 

b.  Inertia 

c.  Truck  unbalance 
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SECTION  III 


STOPPING  DISTANCE  DERIVATION 


1.  AIR  AND  TRANSITION  DISTANCE 

Approach,  flare,  and  transition  distances  are  important  when  considering  total  stopping 
distance.  The  major  factors  that  affect  the  actual  distances  involved  in  these  segments  are: 

• Approach  velocity 

• Glide  slope 

• Height  above  threshold 

• Flare  load  factor 

• Time  to  brake  application 

• Aircraft  deceleration 

A digital  program  was  used  to  analyze  the  effect  of  these  parameters  on  stopping  distance. 
Figures  4 through  7 summarize  the  typical  variations  in  distance  that  can  be  expected  from 
a parameter  change.  The  results  sho'-'n  are  general  and  applicable  to  any  aircraft  by 
appropriate  variable  selection.  Figures  4 and  5 point  out  that  glide  slope  and  height  above 
the  threshold  can  cause  the  approach  distance  to  vary  by  as  much  as  250%,  while  the  flare 
load  factor  and  approach  velocity  can  cause  flare  distance  variations  or  300%.  The  transition 
distance  is  primarily  a function  of  approach  speed,  aircraft  deceleration,  and  time  to  brake 
application.  Figure  6 indicates  that  delaying  brake  application  can  result  in  the  loss  of 
significant  runway  length,  however  reasonable  variations  in  approach  speeds  do  not  result  in 
large  distance  changes.  Figure  7 shows  that  reasonable  changes  in  aircraft  deceleration  do 
not  cause  significant  distance  variations. 

The  six  factors  previously  listed  are  only  a partial  list  ot  parameters  affecting  the  air  and 
transition  distances.  Other  influencing  parameters  involve  aerodynamic  values,  the  landing 
gear  system,  airplane  configuration,  and  environmental  factors.  These  parameters  do  not, 
however,  cause  significant  variations  in  air  and  transition  distances.  The  significant  factors 
(approach  velocity,  glide  slope,  height  above  the  threshold,  flare  load  factor,  and  time  to 
brake  application)  are  parameters  that  reflect  pilot  technique.  Because  pilot  technique  is  the 
major  factor  involved  in  the  determination  of  air  and  transition  distance,  the  prediction  of 
total  stopping  distance  is  difficult;  therefore,  the  only  parametric  changes  pursued  during 
this  program  were  those  involved  in  the  braking  portion  of  the  stop.  The  data  scatter 
resulting  from  variations  in  pilot  technique  can,  however,  be  used  to  develop  advisory 
information  for  aircraft  operation,  which  was  outside  the  scope  of  this  program. 


11 


(ft)  Approach  distance  (ft) 


Figure  ^.-Approach 
Distance  us 
Glide  Slope 


TA  = time  to  brake  application  (sec) 
a = aircraft  deceleration  (ft/sac2) 


Figure  6.~Transition  Distance 
us  Approach  Speed  us 
Time  to  Brake  Application 


Figure  S. —Flare  Distance  us 
Approach  Speed 


175  200  225  250 

Approach  velocity  (fps) 


Figure  /.-Transition  Distance 
vs  Approach  Speed  us 
Aircraft  Deceleration 


A summary  of  the  equations  used  to  calculate  the  approach,  flare,  and  transition  distances 
is  given  in  Appendix  A. 


2.  BRAKING  DISTANCE 

The  brake  control  simulation  (Ref.  1)  used  to  calculate  braking  distance  represents  a refine- 
ment of  an  earlier  baseline  antiskid  simulator  used  at  Boeing. 

The  development  of  the  advanced  simulation  is  the  result  of  many  flight  test  programs 
Boeing  has  conducted  in  recent  years.  In  addition  NASA,  FAA,  and  USAF  data  from  recent 
Boeing  727,  F-4,  and  C-141  aircraft  tests  was  reviewed.  Flight  test  data  has  been  directly 
compared  with  simulator  and  airplane  control  system  tests.  The  outgrowth  of  the  program 
has  led  to  many  refinements  and  the  development  of  more  sophisticated  models.  As  a result, 
the  simulation  more  closely  represents  the  airplane  and  can  accurately  evaluate  the  effect  of 
various  aircraft,  brake,  and  skid  control  system  parameters  on  stopping  performance  and 
landing  gear  stability.  The  payoffs  from  this  undertaking  are  already  incorporated  in  the  737 
Advanced,  Advanced  727,  and  Improved  747F  skid  control  systems. 

The  simulator  is  ai.  analog-hardware  system  consisting  of  three  analog  computers  simulating 
the  aircraft  dynamics  and  a brake  hydraulic  system  mockup.  All  actual  hydraulic  system 
parts  that  influence  brake  system  performance  are  incorporated.  ! ■ addition,  a skid  control 
card  or  box  is  used.  In  this  way,  the  electronic  and  hydraulic  characteristics,  which  may 
include  nonlinearities,  are  accurately  reproduced. 

The  computer  simulation  consists  of  five  interrelated  elements;  (1)  airplane  dynamics,  (2) 
strut  dynamics,  (3)  wheel  dynamics,  (4)  brake  torque,  and  (5)  ground  force.  The  inter- 
relation of  these  elements  and  the  hardware  portion  of  the  simululion  are  shown  in  Figure  8. 
The  basic  airplane  dynamics  include  airplane  pitching,  aerodynamic  lift  and  drag,  engine  idle 
thrust,  and  tire-to-ground  forces.  Also  included  are  detailed  simulations  of  the  landing  gear 
strut  and  the  tire-wheel  dynamics.  Multiple  wheel  brukitig  is  also  provided  in  the  simulation. 

Although  the  actual  brake  is  used  to  determine  the  pressure-volume  characteristics,  the 
pressure  torque  relationship  is  simulated  on  the  computer  using  the  pressure  feedback 
information  from  the  actual  brake,  The  brake  simulation  also  includes  steady-state  torque 
gain,  retractor  spring  deadband,  torque  fade,  torque  peaking,  and  frequency  response.  The 
analytical  models  used  were  derived  as  a result  of  many  hardware  tests  conducted  at  brake 
vendor  facilities  and  tests  conducted  at  AFFDL  Landing  Gear  Test  Facility  on  two  different 
occasions. 

A detailed  description  of  the  Boeing  Brake  Control  Simulation  is  given  in  Appendix  B. 
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SECTION  IV 

SENSITIVITY  STUDY  TEST  PROCEDURE 


The  objective  of  the  sensitivity  study  was  to  determine  the  parameters  that  influence  stop- 
ping distance.  To  properly  assess  the  parameters  and  draw  some  general  conclusions,  five 
aircraft  were  studied.  The  Lockheed-Georgia  C-141A,  McDonnell-Douglas  F-4E,  Boeing 
747-200,  Boeing  737  Advanced,  and  Boeing  727-200  were  chosen,  with  NASA,  FAA,  and 
USAF  concurrence. 

The  data  used  in  the  brake  control  simulator  for  the  Boeing  airplanes  (727,  737,  747)  is  well 
documented  within  The  Boeing  Company.  It  represents  a compilation  of  the  currently 
accepted  data  used  in  various  Boeing  airplane  simulators  and  was  checked  for  accuracy  to 
ensure  the  reliability  of  the  results  obtained  from  the  computer  study.  Data  for  the  F-4  and 
the  C-141  were  obtained  from  the  USAF.  This  data  was  reduced  to  the  form  required  for 
use  in  the  simulation.  Note  that  this  report  generally  uses  a shortened  form  of  the  aircraft 
designations.  As  used,  the  short  terms  are  synonymous  with  the  specific  model  designations. 

As  a starting  point,  a baseline  airplane  was  defined  for  each  aircraft.  The  baseline  airplane 
represents  an  aircraft  in  a three-point  taxi  attitude  and  of  typical  (mid-range)  landing 
weight,  approach  speed,  center-of-gravity  location,  landing  flap  setting,  and  engine  thrust. 
The  actual  parameters  required  for  the  airplane  simulator  are  defined  in  Table  2.  Table  3 
lists  the  baseline  values  used  for  each  airplane. 

During  the  sensitivity  study  each  parameter  was  changed  and  the  new  value  of  braking 
distance  was  evaluated  on  the  simulator.  The  range  over  which  a parameter  was  varied 
reflected  values  observed  in  normal  service  of  the  airplane.  Some  of  the  variables  were  not 
independent,  and  those  groups  of  interrelated  parameters  were  varied  appropriately  to- 
gether. An  example  of  this  is  stall  speed  and  gross  weight. 

The  results  of  the  parametric  study  »»  lid  for  the  purpose  of  this  program  and  contribute 
a whole  new  dimension  to  the  state  of  knowledge  for  developing  more  efficient  stopping 
systems  and  more  precise  methods  of  predicting  aircraft  stopping  performance.  The 
parameters  are,  however,  not  varied  in  this  manner  in  real  life;  therefore,  the  results  should 
not  be  misused  to  judge  operational  performance  of  any  aircraft  studied. 


IS 


Table  2. -Definition  of  Simulator  Parameters 


Para  mater 

Definition 

Unit 

Airplane  Parameter! 

. AW 

Effective  wing  area 

ft* 

CD 

Drag  coefficient 

COD 

Drag  coefficient  without  spoilers 

CL 

Lift  coefficient 

CLD 

Lift  coefficient  without  spoilers 

FEO 

Engine  Idle  Thrust  at  zero  velocity 

Ibf 

HB 

Height  of  CG  above  ground 

ft 

IYV 

Muss  moment  of  Inertia,  pitch 

ft-lb-sac^ 

KE 

Change  of  Idle  thrust  with  velocity 

lbf'sec/ft 

LA 

Nose  gear  to  CG  distance 

" 

LB 

Main  gear  to  CG  distance 

ft 

NB 

Number  of  brakes  per  main  strut 

NBA 

Number  of  moin  gear  brakes  per  airplane 

NBN 

Number  of  nose  gear  wheels 

NS 

Number  of  m«ln  dear  ftruts  par  airplane 

RHO 

Air  density 

ibf-sec^/ft4 

VI 

Initial  airplane  velocity 

f ./«c 

VSTOP 

Final  airplane  velocity 

ft/sec 

WA 

Weight  of  airplane 

Ibf 

Brake  Parameter! 

KP 

Torque  peeking  gain 

MB 

Mass  of  brake  heat  sink 

Ibm 

OMGP 

Wheel  velocity  at  start  of  torque  peaking 

rad /sec 

PC 

Retractor  spring  pressure 

psi 

TBG 

Torque  gain 

ft-lbf/psi 

THBB 

Temperature  at  initiation  of  fade 

°F 

WN 

Natural  frequency  of  torque  respond 

Hz 

ZETA 

Damping  ratio  of  torque  response 

Tire  Parameter* 

D 

Tire  diameter 

in. 

DO 

Tire  deflection 

in. 

IW 

Mass  moment  of  inertia  of  tire,  wheel,  and  brake 

ft-lbf-sec^ 
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Table  2. -Definition  of  Simulator  Parameters  (Concluded) 


Parameter 

Definition 

Unit 

Tlra  Pirtmtttfi  (Continued) 

PI 

Tlra  oparatlng  Inflation  prauure 

psi 

PR 

Tire  rated  inflation  pressure 

Pii 

RR 

Tire  rolling  radius 

ft 

RT 

Tire  torque  radius 

ft 

Strut  Paramatart 

CO 

Main  gear  vertical  damping  coefficient 

Ibf -tec/ft 

CON 

Nose  gear  vertical  damping  coefficient 

Ibf-Mc/ft 

cs 

Main  gear  fore-aft  damping  coefficient 

lbf-sec/ft 

CT 

Torsional  Damping  between  strut  and  brake 

Ibf-ft-sec 

IS 

Mass  moment  of  Inertial  of  main  gear  strut 

ft-lbf-sec^ 

KO 

Main  gear  vertical  stiffness 

lbf/ft 

KON 

Nose  gear  vertical  stiffness 

Ibf/ft 

KS 

Main  gear  fore-aft  stiffness 

lbf/ft 

L 

Effective  strut  length 

ft 

MS 

Effective  mass  of  strut 

Ibf-sec^/ft 

Table  3. -Baseline  Parameter  Values  Used  in  Airplane  Simulation 
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SECTION  V 

SENSI1  VITY  CTUDY  TEST  CONDITIONS 


The  sensitivity  study  involved  changing  a parameter  or  group  of  parameters  to  analyze  its 
influence  on  stopping  distance.  Based  on  the  factors  listed  in  Table  1 , a set  of  test  condi- 
tions was  formulated.  The  test  conditions,  as  listed  In  Table  4,  were  performed  on  the 
analog-hardware  simulator  to  analyze  braking  distance.  Table  S lists  the  parameters  changed 
and  the  numerical  value  of  the  parameter  associated  with  each  test  condition  for  the  five 
aircraft. 

A test  outline  was  proposed  to  NASA,  FAA,  and  USAF  and,  with  their  concurrence,  nine 
sensitivity  tests  were  formulated  to  analyze  the  sensitivity  of  the  various  aircraft  systems  to 
a parameter  change.  The  tests  were  performed  selectively  at  each  test  condition.  The  nine 
sensitivity  tests  have  been  divided  into  three  major  categories.  The  categories  and  the  tests 
associated  with  each  are  as  follows: 

• Stability  studies: 

Test  I -strut  stability 

• Performance  studies: 

Test  2-touchdown  dynamics 
Test  3 -stabilized  landing 
Test  4-mu  steps 
Test  5 -wet  runway 

• Hydraulic  system  studies: 

Test  6-frequency  response 
Test  7-step  response 

Test  8 -antiskid  valve  characteristics 

Test  9-brake  pressure  and  volume  characteristics 

A detailed  description  of  the  test  procedure  and  sequence  can  be  found  in  ASD-TR-74-41 , 
Volume  II,  Section  IX.  The  three  majot  categories  of  sensitivity  tests  are  briefly  described 
below  to  point  out  their  general  significance. 
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4. -Test  Conditions 

Baseline  Study 

Nominal  values  of  all  parameters 
Parametric  Studies 
Airplane 


1 . Weight 

a.  Maximum  landing 

b.  Minimum  landing 

2.  Center  of  Gravity 

a.  High 

b.  Low 

c.  Forward 

d.  Aft 

3.  Brake  application  speed 

a.  +10% 

b.  +20% 

4.  Aerodynamics 

a.  Spoiler  or  drag  chute  deployment  - 1 .0  sec  after  touchdown 

b.  Spoiler  or  drag  chute  deployment  - 2.0  sec  after  touchdown 

c.  No  spoilers  or  drag  devices 

d.  60%  effective  spoilers 

e.  40%  effective  spoilers 

f.  120%  engine  idle  thrust 

g.  80%  engine  idle  thrust 

6.  Pilot  technique 

a.  150%  nominal  pressure  application  rate 

b.  60%  nominal  pressure  application  rate 

c.  Nominal  rate  at  2.0  sec  from  touchdown 

d.  Nominal  rate  at  4.0  sec  from  touchdown 

e.  75%  of  full  metered  pressure 

f . 60%  of  full  metered  pressure 

Runwa*  <nd  environmental  system 

1.  Wind 

a.  10  knots 

b.  20  knots 

c.  -10  knots 


Table  4,-Test  Conditions  (Continued) 


Parametric  Studies  (Continued) 

Runway  and  environmental  system  (continued) 

2. 

Air  density 

s. 

Hot  day  (83°F  - 28°C) , High  altitude  (6000  ftl 

b. 

Cold  day  <-60°F  - -61°CI  tea  level 

3. 

Runway 

a. 

Roughness  profile 

Landing  gear  systems 

1. 

Brake 

torque  characteristics 

a. 

High  fade  brake  - 80%  MB  (only  0.4  mu  and  greater) 

b. 

Low  fade  brake  - 00%  MB  (only  0.4  mu  and  greater) 

c. 

Peaking  to  1 50%  of  running 

d. 

No  peaking 

e. 

Torque  response  breakpoint  - 160%  of  nominal 

f. 

Torquj  response  breakpoint  - 50%  of  nominal 

9. 

Torque  gain  120%  of  nominal 

h. 

Torque  gain  80%  of  nominal 

1. 

Variable  torque  gain  T - f(p)®-5 

I. 

Linear  torque  gain  T e f(p) 

2. 

Tire 

a. 

Inflation  pressure  120%  of  nominal 

b. 

Inflation  pressure  80%  of  nominal 

c. 

50%  worn  tire 

d. 

80%  worn  tire 

e. 

Low  tire  heating 

f. 

Flat  m-  apeak 

3. 

Strut 

a. 

Maximum  strut  frequency  varying  mass 

b. 

Minimum  strut  frequency  varying  mass 

c. 

Maximum  strut  frequency  varying  stiffness 

d. 

Minimum  strut  frequency  varying  stiffness 

B. 

Vertical  stiffness  120%  of  nominal 

f. 

Vertical  stiffness  80%  of  nominal 

0- 

Vertical  damping  120%  of  nominal 

h. 

Vertical  damping  80%  of  nominal 
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Table  4. -Test  Conditions  (Concluded) 

Parametric  studies  (continued) 

Hydraulic  syttem 

a.  Decrease  line  diameter  50%  (nearest  diameter) 

b.  increase  line  diameter  50%  (nearest  diameter) 

c.  Move  dynamic  breakpoint  out  1 50%  of  nominal  (line  length  change) 

d.  Move  dynamic  breakpoint  in  60%  of  nominal  ( line  length  change) 

e.  Insert  20%  restriction  in  return  line 

f.  Increase  brake  volume  by  10  cu  in. 

g.  Increase  brake  p-v  gain 


m 


Table  5. -Parameter  Change  values  (Continued) 


Teat  condition  and 
parameter  changed 


Airplane  model 


Landing  Gear  Syttemi 
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Values  (Concluded) 


1.  STABILITY  STUDIES 

System  stability  is  directly  related  to  stopping  performance.  Severe  instability  can  result  in 
the  loss  of  braking  and  can  cause  serious  safety  hazards.  The  study’s  purpose  was  to  evaluate 
the  ability  of  a brake  control  system  to  contribute  to  the  stability  of  the  gear . 

2.  PERFORMANCE  STUDIES 

The  performance  studies  provide  a measure  of  the  performance  capability  of  the  brake 
system.  The  tests  performed  fall  into  two  categories.  The  first  defines  the  operation  of  the 
airplane  under  stable  landing  conditions.  The  second  evaluates  the  ability  of  the  brake 
system  to  adapt  to  the  typical  dynamic  operating  conditions  encountered  during  an  actual 
landing. 

3.  HYDRAULIC  SYSTEM  STUDIES 

The  hydraulic  system  studies  measure  the  response  of  the  antiskid  valve,  control  box,  and 
the  actual  brake  hydraulic  system.  Specific  tests  were  designed  to  define  both  the  overall 
and  component  performance  of  the  system.  The  results  provide  an  insight  into  aircraft 
braking  system  performance  and  can  be  used  to  further  improve  some  of  the  systems. 


SECTION  VI 
TEST  RESULTS 


I.  BASELINE 

A baseline  airplane  was  defined  for  each  of  the  five  aircraft  studied.  The  baseline 
configuration  is  meant  to  represent  a typical  aircraft  during  its  landing  phase.  The  numerical 
values  of  the  baseline  parameters  are  listed  in  Table  3.  The  values,  when  used  in  the  brake 
control  simulation,  establish  a unique  distance  versus  ground  friction  relation  for  each 
airplane. 

Figure  9 relates  the  braking  distance  and  peak  available  ground  friction  (peak  available  mu) 
for  each  of  the  five  baseline  airplanes.  The  data  is  associated  with  the  braking  segment  only, 
so  the  distances  shown  represent  only  braked  distances;  'pproach,  Hare,  and  transition 
distances  are  excluded.  For  a skidding  tire,  the  instantaneous  coefficient  of  friction  obtained 
at  the  tire-runway  Interface  depends  on  the  condition  of  the  runwuy,  aircraft  characteristics, 
tire  properties,  and  tire  slippage.  Thus,  peak  available  mu  (such  as  shown  in  Figure  9)  is  a 
computer  Input  that  defines  the  muxlmum  value  of  friction  available  between  the  tire  and 
ground  during  a test  condition. 

An  ideal  brake  system  should  operate  at  the  peak  available  friction  value  during  the  entire 
stop.  However,  the  actual  anMsxld  efficiencies  realized  are  the  result  of  component  sizing 
und  system  characteristics  and  Integration.  In  addition,  some  compromises  are  necessary,  so 
the  system  actually  functions  over  a range  of  fr-tlon  somewhat  lower  than  the  peak  value. 
Since  the  actual  antiskid  system  has  been  used  in  the  brake  control  simulation  the  distances 
produced  (Figure  9)  reflect  the  efficiency  of  the  braking  system. 

2.  SENSITIVITY 

The  sensitivity  study  of  the  five  airplanes  involved  changing  a parameter,  or  group  of 
parameters,  and  observing  the  effect  on  braking  distance.  The  bruking  distance  associated 
with  a change  was  then  compared  to  the  baseline  airplane  distance  at  the  same  value  of  peak 
available  mu.  In  this  manner,  the  effect  of  a specific  change  could  be  analyzed 
quantitatively. 

The  braking  data  obtained  from  the  brake  control  simulator  represents  an  absolute  distance. 
To  facilitate  the  analysis  of  a parameter  change,  a normalized  distance  has  been  introduced, 
it  Is  termed  "baseline  braking  distance  percentage.”  The  use  of  a normalized  distance  allows 
the  five  airplanes  to  be  conrpured  simultaneously.  The  baseline  braking  distance  percentage 
Is  the  bruking  dUance  associated  with  a parameter  change  divided  by  the  baseline  braking 
distance  at  the  same  value  of  peak  available  mu;  this  quotient  is  then  multiplied  by  100%  to 
obtain  a percentage  value. 

Each  airplane  uses  its  own  baseline  distances  as  normalizing  factors.  Distances  longer  than 
baseline  are  greater  than  1 00%,  and  shorter  distances  are  reflected  as  less  than  1 00%. 


The  braking  distance  results  from  the  sensitivity  study  have  been  reduced  to  bar  charts, 
Figure  10,  Sheets  1 through  15.  The  test  condition  corresponding  to  each  bar  chart  is  given 
to  the  left  and  above  the  chart.  On  each  graph,  the  baseline  braking  distance  percentage  for 
the  five  airplanes  has  been  plotted  for  three  values  of  peak  available  mu:  0.6, 0.4,  and  0.2. 

A brief  analysis  of  the  results  follows.  The  raw  data  along  with  the  associated  performance 
indices  can  be  found  in  ASD-TR-74-41 , Volume  II  Section  X. 

a.  AIRPLANE  FLIGHT  CHARACTERISTICS 

1.  Landing  Weight  (la  and  lb,  Sheet  1,  Figure  10) 

Three  parameters  were  changed  in  these  tests:  landing  weight,  landing  speed,  and  mass 
pitching  moment.  The  landing  speed  was  adjusted  to  reflect  changes  in  stall  speed  resulting 
from  variations  in  landing  weight.  The  mass  pitching  moment  was  varied  to  reflect  changes 
in  load  distribution.  As  shown  in  the  bar  charts,  braking  distance  increases  with  an  increase 
in  gross  weight.  This  results  from  the  increased  kinetic  enetgy  that  must  be  dissipated  during 
braking.  The  variation  in  distance  is  largely  attributed  to  the  increase  in  brake  application 
speed.  Additional  tests  indicated  that  the  weight  variation  alone  causes  about  a 3%  variation 
in  the  baseline  braking  distance. 

2.  Center  of  Gravity 

(a)  Htyh  and  Low  Center  of  Gravity  (2a  and  2b,  not  diagramed).-  High  and  low  center  of 
gravity  tests  were  run  on  the  Boeing  737.  Initial  digital  computer  studies  showed  this  factor 
to  be  insignificant  to  braking  distance,  and  actual  tests  with  the  737  aircraft  verified  the 
prediction.  As  a result,  these  tests  were  deleted  for  the  four  remaining  airplanes. 

(bl  Forward  and  Aft  Center  of  Gravity  (2c  and  2d,  Sheet  1,  Figure  10). -Results  show  that 
aft  location  of  center  of  gravity  (CG)  causes  a decrease  in  braking  distance.  This  results 
because  more  weight  is  placed  on  the  main  gear,  increasing  the  available  braking  force.  When 
analyzing  CG  location,  steering  requirements  must  also  be  considered.  An  increase  in  main 
gear  load  will  result  in  loss  of  rose  gear  load  and  thus  cornering  capability.  On  wet  runways, 
this  situation  warrants  serious  consideration,  and  some  compromise  is  required. 

The  F-4  is  the  notable  exception  to  the  general  trend.  This  airplane  uses  the  Hydro-Aire 
Mark  II  system  similar  to  the  Boeing  727,  though  the  choice  of  thresholds  in  the  F-4 
antiskid  circuitry  is  different.  As  a consequence  of  this  and  system  tuning,  the  F-4  antiskid 
system  operates  predominatly  on  the  front  side  of  the  mu-slip  curve.  With  more  weight  on 
the  main  gear,  increased  cycling  of  brake  pressure  occurs,  which  decreases  system  efficiency 
(mean  effective  pressure)  and  results  in  increased  stopping  distances.  On  the  other  hand,  a 
forward  CG  location  resulting  in  less  weight  on  the  main  gear  causes  the  system  to  cycle  less 
often,  thus  increasing  the  system  efficiency. 

3.  Brake  Application  Speed  (3a  and  3b,  Sheet  2,  Figure  10) 

The  brake  application  speed  was  varied  to  show  the  effect  of  a pilot  making  a “hot"  landing. 
The  results  reflect  the  additional  braking  distance  required  to  dissipate  the  increased  kinetic 
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Figure  W.-Braking  Distance  Test  Results  (Sheet  7 of  15) 
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Figure  10,-Braking  Distance  Test  Results  (Sheet  8 of  15i 
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c.  Da  create  Line  Length 


Figure  10. -Braking  Dis 
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Test  Results  (Sheet  14  of  15) 


energy  of  the  aircraft.  The  F-4  appears  to  operate  quite  differently  at  low  ground  friction 
values,  largely  because  of  its  lift  and  drag  characteristics.  An  increase  in  the  landing  speed  of 
the  F-4  results  in  a large  increase  in  drag  and  a small  decrease  in  main  gear  loading.  At  the 
low  mu  condition,  the  increased  kinetic  energy  is  nearly  compensated  by  the  increase  in 
drag.  Thus  the  braking  distance  shows  only  slight  variations. 

4.  Spoiler  or  Drag  Chute  Deployment  (4a  and  4b,  Sheet  2,  Figure  1 0) 

The  general  trend  shown  is  that  delaying  spoiler  application  or  drag  chute  deployment,  is 
deterimental  to  braking  performance.  Significant  braking  force  is  lost  during  the  period 
prior  to  spoiler  deployment.  From  a braking  performance  standpoint,  the  pilot  should  spoil 
the  lift  as  soon  as  possible. 

5.  Spoiler  or  Drag  Chute  Effectiveness  (4c,  4d,  and  4e,  Sheet  3,  Figure  10) 

As  shown,  braking  distance  increases  as  the  spoiler  becomes  less  efficient.  The  increased 
distance  results  from:  (1)  loss  of  braking  force  due  to  an  increase  in  life,  and  (2)  loss  of 
effective  drag  force  due  to  lower  drag  coefficients.  The  graphs  show  the  increasing  effect  of 
spoiler  effectiveness  at  low  mu  conditions.  This  reflects  the  importance  of  spoilers  during 
adverse  weather. 

6.  Engine  Idle  Thrust  (4f  and  4g,  Sheets  3 and  4,  Figure  1 0) 

Changes  in  the  engine  idle  thrust  value  generally  cause  a change  in  braking  distance  of  less 
than  3%.  The  F-4  is  the  exception;  it  has  an  unusually  high  ratio  of  idle  thrust  to  weight. 
Thus  a change  in  the  idle  thrust  value  causes  a significant  change  in  the  kinetic  energy  to  be 
dissipated  during  braking.  This  in  turn  directly  affects  the  braking  distance. 

7.  Brake  Pressure  Application  Rate  (5a  and  5b,  Sheet  4,  Figure  10) 

Distance  variations  resulting  from  a change  in  pressure  application  rate  are  a maximum  of 
4%  and  in  general  average  about  1%.  For  the  range  of  pressure  application  rates  tested,  it 
does  not  appear  that  the  parameter  has  much  effect. 

8.  Delay  of  Brake  Pressure  Application  (5c  and  5d,  Sheets  4 and  5,  Figure  10) 

The  delayed  pressure  application  bar  charts  include  both  the  actual  braking  distance  and  the 
added  transition  distance  caused  by  the  delay,  As  can  be  seen,  the  effect  of  delaying  brakes 
simply  causes  an  increase  in  total  ground  distance.  This  results  from  the  loss  of  braking 
during  the  delay  period  and  the  small  amount  of  speed  lost  during  that  time. 

9.  Metered  Pressure  (5e  and  5f,  Sheet  5,  Figure  10) 

The  results  from  the  metered  pressure  tests  do  not  show  a general  distance  trend.  However, 
the  variation  of  braking  distance  between  the  different  aircraft  is  attributable  to  valve 
characteristics,  torque  limiting  of  the  brake,  and  skid  control  adaptation  to  conditions.  The 
Boeing  747,  737,  and  727  show  an  increase  in  braking  distance  as  the  metered  pressure  is 
lowered.  For  these  aircraft,  the  brake  is  pressure-limited  and  cannot  develop  sufficient 
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torque  to  initiate  a skidding.  The  F-4  and  C-141  show  an  opposite  trend.  At  lower  pressures, 
the  dynamic  response  of  these  systems  is  slower,  reducing  the  cycling  of  the  antiskid  system, 
increasing  system  efficiency,  and  reducing  stopping  distances. 


b.  RUNWAY  AND  ENVIRONMENTAL  SYSTEM 

1.  Wind  (la,  lb,  and  lc.  Sheet  6,  Figure  10) 

The  environmental  parameter  wind  has  quite  a straightforward  effect  on  braking  distance. 
Headwinds  (positive-knot  values)  reduce  braking  distance;  tailwinds  increase  the  distance. 


2.  Hot  and  Cold  Daysdd  and  le,  Sheet  7,  Figure  10) 

During  these  tests,  the  air  density  was  varied  to  change  the  lift  and  drag  forces  on  the 
airplane.  On  hot  days,  both  lift  and  drag  decrease;  on  cold  days,  they  increase.  The  results 
show  generally  that  drag  has  the  greater  effect  on  braking  distance.  On  hot  days,  the 
distances  are  longer  because  of  the  decreased  aerodynamic  drag.  The  Boeing  747  is  an 
exception;  its  lift-to-drag  ratio  is  much  greater  than  the  other  four  airplanes,  and  the  lift 
effect  is  more  significant.  During  a hot  day,  lift  decreases,  placing  more  weight  on  the  main 
gear.  This  increases  the  braking  capability  and  decreases  the  braking  distance.  In  operational 
situations,  higher  flap  settings  are  selected  to  compenstate  for  this  effect. 

3.  Rough  Surface  (3a,  Sheet  7,  Figure  10) 

The  addition  of  the  rough  runway  surface  has  the  effect  of  providing  an  additional  perturba- 
tion to  the  vertical  motion  of  the  aircraft.  This  results  in  random  main  gear  loading  and  tests 
the  capability  of  a system  to  adapt  to  adverse  runway  surfaces.  The  data  indicates  that  each 
system  tested  had  no  problem  adapting  to  this  surface  condition.  The  power  spectral  density 
curve  of  the  rough  runway  condition  simulated  is  given  in  Appendix  B. 

c.  LANDING  GEAR  SYSTEMS 

I.  Brake  FadeOa  and  lb,  Sheet  8,  Figure  10) 

During  the  high  and  low  brake  fade  tests,  the  effective  mass  of  the  brake  was  lowered  to 
initiate  brake  fade  earlier  in  the  braking  segment.  The  largest  variation  in  braking  distance 
from  this  parameter  change  is  2%.  In  light  of  the  Idf  repeatability  accuracy  of  the  simulator, 
brake  fade  is  not  considered  a significant  variable  in  braking  performance  during  landing 
rollout. 
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2.  Torque  Peaking  (lc  and  id.  Sheet  8,  Figure  10) 

The  nominal  value  of  torque  during  peaking  was  increased  to  1 50%  of  running  torque  and 
decreased  to  no  torque  peaking.  Peaking  is  a low-speed  phenomenon  that  occurs  during  a 
deep  skid  and  at  the  end  of  the  braking  segment.  As  a result,  torque  peaking  has  little  or  no 
effect  on  the  actual  braking  distance. 

3.  Torque  Response  Frequency  (le  and  If,  Sheet  8,  Figure  10) 

A change  in  the  torque  response  frequency  breakpoint  effects  the  lag  in  the  pressure-torque 
model.  As  predicted  by  theory  and  borne  out  by  the  data,  movement  of  the  breakpoint  out 
to  150%  of  nominal  reduces  the  system  lag  and  allows  faster  response  and  shorter  braking 
distances.  Conversely,  movement  of  the  breakpoint  inward  increases  the  lag  and  lengthens 
the  distance.  This  aspect  of  brake  performance  should  be  stressed  to  brake  suppliers. 

4.  Torque  Gain  (Ig,  ih,  li,  and  U,  Sheet  9,  Figure  10) 

During  these  four  tests,  the  brake  pressure-torque  relationship  was  changed.  The  airplanes 
show  a random  degree  of  sensitivity  to  the  changes.  Little  data  is  available  on  the  actual 
pressure-torque  relationship  of  brakes.  The  observed  sensitivity  simply  points  to  the  need 
for  a better  understanding  of  the  torque  gain  of  a brake. 


5.  Mu-Slip  Curve  Variations  (2a,  2b,  2e,  and  2f,  Sheet  10,  Figure  10) 


Variations  in  the  shape  of  the  mu-slip  curves  were  made  to  analyte  how  each  antiskid 
system  would  react  to  changes  in  runway  traction.  The  data  shows  significant  sensitivity, 
which,  however,  is  random.  The  randomness  between  the  different  aircraft  results  from  the 
manner  in  which  the  antiskid  system  is  implemented.  Each  system  and  the  mu-slip  curve 
variation  f analyzed  in  depth  in  ASD-TR-7‘‘41 , Volume  II,  Section  VIII.  The  general 
conclusion  to  be  made  from  these  data  is  that  the  mu-siip  curve  shape  influences  stopping 
distance. 


6.  Worn  Tire  (2c  and  2d,  Sheet  1 1,  Figure  10) 

During  the  worn  tire  tests,  the  rolling  radius,  torque  radius  and  wheel-tire  inertia  were 
changed.  The  change  in  braking  distance  averaged  about  1%  exclusive  of  the  F-4.  The  F-4  is 
unusual  because  of  its  small  main  gear  tire.  The  fixed  change  in  torque  r-'-Mus  (used  for  all 
airplanes)  represents  a large  reduction  in  ground  torque.  This  reduced  to.  -]ue  has  caused  a 
significant  increase  in  braking  distance, 

7.  Strut  Frequency  (3a,  3b,  3c,  and  3d,  Sheet  12,  Figure  10) 

The  natural  frequency  of  the  main  gear  was  changed  to  check  for  the  existance  of  coupling 
between  the  skid  control  box  logic  and  gear  motion.  When  coupling  exists,  a catastrophic 
gear  failure  may  occur,  or  braking  performance  may  be  jeopardized.  The  727  is  the  only 
aircraft  where  the  braking  distance  was  effected  by  greater  than  2%.  In  this  case,  the 
coupling  lowered  the  mean  antiskid  pressure  level.  This  reduced  the  average  brake  torque 
throughout  the  run,  which  increased  the  braking  distance. 
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8.  Vertical  Stiffness  and  Vertical  Damping  (3e,  3f,  3g,  and  3h,  Sheet  12,  Figure  1C 

During  these  tests,  the  characteristics  of  the  main  and  nose  gear  were  changed.  Such  a 
change  primarily  effects  the  pitching  motion  of  the  aircraft  and,  to  some  extent,  the  tire 
load.  The  data  shows  braking  distance  to  be  insensitive  to  vertical  gear  stiffness  and  damping 
variations. 

d.  HYDRAULIC  SYSTEMS 

To  realize  the  best  possible  stopping  performance,  each  subsystem  must  be  integrated  into 
the  total  system.  The  hydraulic  system  is  one  such  system  that  can  be  tuned  to  aid  in 
obtaining  optimum  performance.  These  tests  were  designed  to  indicate  how  well  the  hydrau- 
lic system  had  been  optimized.  No  attempt  has  been  made  to  retune  the  antiskid  system 
during  these  tests.  Decreased  braking  performance  may  be  reclaimed  by  retuning  the  anti- 
skid control  box. 

1.  Line  Diameter  (a  and  b,  Sheet  14,  figure  10) 

Line  diameter  tests  reveal  that,  with  the  exception  of  the  F-4,  most  of  the  systems  have 
been  properly  implemented.  The  F-4  data  shows  that  substantial  performance  can  be  gained 
by  increasing  the  hydraulic  line  size.  Decreased  line  diameter  tests  were  not  performed  on 
the  F-4  and  747  because  the  standard  lines  were  I /4-in.  diameter,  and  further  reduction  was 
not  deemed  feasible. 

2.  Line  Length  (c  and  d.  Sheet  1 4,  Figure  1 0) 

Changes  in  line  length  effect  the  frequency  response  gain  and  lag  of  the  hydraulic  system. 
Basic  hydraulic  analysis  predicts  that  shorter  lines  decrease  the  gain  and  lug  in  the  system; 
this  will  result  in  better  performance.  Excluding  the  F-4,  the  frequency  response  and  brak- 
ing distance  showed  the  predicted  trend.  The  F-4,  however,  exhibited  a decrease  in  gain  as 
line  length  increases,  thus  braking  distance  decreases  as  length  increases.  Further  turning  of 
the  skid  control  system  or  use  of  a different  skid  control  system  and  associated  hardware 
could  have  shown  different  results.  This,  however,  was  not  within  the  scope  of  these  studies. 

3.  Restriction  in  Return  Line  (e.  Sheet  15,  Figure  10) 

Except  for  isolated  cases,  braking  distance  was  unaffected  by  hydraulic  return  line 
restructions. 

4.  Brake  Volume  Changes  (f  and  g,  Sheet  1 5,  Figure  10) 

The  volume  of  the  standard  aircraft  brake  was  changed,  resulting  in  a modification  of  the 
pressure-volume  curve  of  the  brake.  Increasing  the  brake  volume  decreases  the  brake  gain. 
To  increase  the  brake  gain,  the  actual  brake  was  replaced  with  a small  accumulator.  A 
varying  degree  of  sensitivity  was  exhibited  during  these  tests;  however  in  general,  only  the 
increased  gain  case  caused  significant  change.  In  the  case  of  the  747,  the  change  actually 
caused  the  system  to  be  unstable,  and  the  test  could  not  be  performed.  Early  simulator 
development  is  considered  necessary  to  ensure  proper  brake  characteristics.  These  character- 
istics can  be  incorporated  in  brake  specifications.  This  procedure  is  routinely  done  during 
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the  development  of  Boeing  airplanes.  The  F-4  antiskid  system  was  unable  to  accurately 
control  the  brake  pressure.  The  result  was  increased  cycling  of  the  system  and  lowered 
performance.  The  instability  of  the  747  and  the  lower  performance  of  the  F-4  was  caused 
by  an  increase  in  antiskid  system  gain  (resulting  from  the  brake  change)  without  any  tuning 
of  the  antiskid  circuit.  Due  to  proper  tuning  and  hardware  selection  no  stability  problems 
exist  in  the  Boeing  747. 

e.  COMPOSITE  VARIATION  OF  PARAMETERS 

During  the  braking  sensitivity  study,  the  baseline  value  of  a parameter  or  group  of  related 
parameters  was  changed.  The  braking  distance  resulting  from  a change  has  been  used  as  a 
measure  of  that  parameter's  effect  on  braking.  Because  nonlinear  dynamic  equations  are 
involved,  the  effect  of  making  composite  parameter  changes  cannot  be  predicted  by  simple 
algebraic  addition  of  the  effect  of  each  individual  parameter.  This  aspect  is  stressed  to 
discourage  the  use  of  data  presented  for  purposes  other  than  intended  for  this  study. 


SECTION  VII 

PARAMETER  EVALUATION  CRITERIA 

1.  PERFORMANCE  INDICES 

To  evaluate  the  performance  of  a system,  various  normalized  indices  were  developed.  The 
use  of  such  indices  aids  in  the  comparison  of  various  systems. 

The  performance  parameters  used  in  evaluating  each  skid  control  system  were: 

• Airplane  braking  distance 

• Braking  distance  efficiency 

• Developed  Mu  efficiency 

• Skid  index 

• Cornering  index 

With  the  exception  of  braking  distance  efficiency,  these  indices  were  calculated  directly  on 
the  simulator  and  measured  with  a digital  voltmeter.  Braking  distance  efficiency  was  calcu- 
lated from  the  results  of  the  simulator  and  a digital  computer  program. 

a.  AIRPLANE  BRAKING  DISTANCE 

Airplane  braking  distance,  X^,  as  measured  on  the  analog  computer  is  the  distance  the 
airplane  travels  from  brake  application  to  a low-velocity  turn-off  speed. 

b.  BRAKING  DISTANCE  EFFICIENCY 

Braking  distance  efficiency,  ps,  is  the  ratio  of  the  perfect  stop  distance  to  the  braked 
airplane  distance. 


X = VXA  * l00% 


where : 

t?s  = braking  distance  efficiency 
Xp  = perfect  braking  distance 
X^  = airplane  braking  distance 


si  Preceding  pege  blank 


The  perfect  braking  distance  is  the  distance  required  to  stop  the  airplane  if  it  is  braked  for 
the  entire  stop  with  maximum  available  braking  force.  This  distance  is  determined  by 
stopping  a reference  airplane  that  is  identical  to  the  braked  airplane. 

Braking  distance  efficiency  Indicates  the  degree  to  which  the  system  meets  its  primary 
requirement  of  stopping  the  aircraft. 

c.  DEVELOPED  MU  EFFICIENCY 

The  developed  mu  efficiency  is  cumulative  measure  of  the  use  of  available  mu.  It  is  deter- 
mined by  dividing  the  time  integral  of  developed  mu  by  the  time  integral  of  available  mu. 
This  parameter  is  useful  in  evaluating  system  performance  but  can  be  misleading  in  that  a 
high  developed  mu  at  high  speed  is  more  effective  than  a high  developed  mu  at  a low  speed. 
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where: 

n0  = developed  mu  efficiency 
t = time 

T = time  to  simulated  stop 
H0  = developed  mu 
= available  mu 


d.  SKID  INDEX 

The  skid  index  is  a measure  of  the  tire  wear  produced  by  the  brake  control  system.  The 
index  is  defined  as  the  integral  of  the  product  of  sliding  velocity  and  ground  force  divided 
and  normalized  by  the  integral  of  elapsed  lime.  This  parameter  is  also  an  indication  of  side 
force  capability. 


SI 
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where: 

SI  = skid  index 
Vg  = sliding  velocity 
Fq  = ground  force 
T = total  time  for  the  stop 


e.  CORNERING  INDEX 

The  cornering  index  is  a measure  of  the  ability  of  the  tire  to  sustain  a side  load.  The 
instantaneous  cornering  index  is  defined  as  a function  of  tire  slip.  The  total  cornering  index 
value  is  the  summation  of  the  instantaneous  values  throughout  the  stop. 


CI=  100 


/T  (I  -o)ndt 
0 


o 


where:  Cl  = cornering  index,  equal  to  zero  for  a wheel  fully  locked 
and  1 00  for  a free-rolling  wheel 

0 = slip  ratio 

n = sensitivity  coefficient  (2.0) 

T = total  time  for  the  stop 

f.  SYSTEM  STABILITY 


In  addition  to  the  five  performance  indices,  a measure  was  made  of  the  ability  of  an  antiskid 
system  to  contribute  to  system  stability.  The  criterion  used  to  determine  system  stability 
was  that  the  system  is  termed  unstable  if  the  main  gear  strut  oscillations  diverge.  The 
stability  of  the  system  was  measured  by  determining  the  damping  ratio  necessary  to  cause 
instability. 


2.  PARAMETER  RATING  SYSTEM 

The  final  step  in  the  sensitivity  analysis  of  aircraft  braking  performance  is  to  rate  each 
parameter.  To  facilitate  the  rating  of  a parameter  change  the  normalized  distance  “baseline 
braking  distance  percentage”  was  used.  The  use  of  this  term  allows  the  five  airplanes  to  be 
analyzed  as  a group.  In  order  to  uniformly  and  quantitatively  rate  a parameter  change,  the 
following  formula  was  used: 
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^ps  mm 


PR1  = 


s 


baseline  braking  distance 
— _ 


percentage)  - 1 00% 


PRI  = parameter  rating  index 
n * total  number  of  data  points  for  3 
particular  parameter  change 
baseline  braking  distance  percentage)  = the 
baseline  braking  distance  percentage 
value  for  the  ith  data  point 


The  parameter  rating  index  (PRI)  as  calculated  above  is  the  average  percentage  deviation 
from  the  baseline  braking  distance.  Thus,  the  value  of  the  PRI  increases  when  a parameter 
change  causes  larger  deviations  from  the  baseline  braking  distance.  The  data  used  in  the 
calculation  of  the  PRI  are  given  in  ASD-TR-74-41,  Volume  II,  Section  XII.  Also  included 
are  the  final  PRI  values  for  the  dry-stabilized  landing  conditions. 
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SECTION  VIII 

SIMULATOR-TO-AIRPLANE  CORRELATION 


Meaningful  results  from  a simulator  can  be  obtained  only  when  correspondence  to  the 
modeled  system  can  be  demonstrated.  The  simulator  used  during  this  study  is  the  culmina- 
tion of  considerable  effort  including  model  development,  subsystem  test  and  evaluation, 
correlation  with  airplane  tests,  and  operational  usage. 

The  analog  computer  models  used  in  the  simulator  are  expanded  versions  of  those  used  in 
the  Boeinc  Baseline  Simulation  (Ref.  1).  A continuous  improvement  effort  has  been  made 
to  represent  the  dynamic  system  behavior  more  thoroughly  as  data  has  become  available. 
Data  used  for  model  expansion  and  improvement  has  been  derived  from: 

• Flight  test  results 

• Dynamometer  tests 

• Contracts  such  as  “Brake  Sensitivity  Analysis”  and  “Antiskid  Performance  Compatabil- 
ity  ind  Improvement” 

• In-house  research 

• Literature 

The  basic  simulation  has  been  used  for  system  evaluation,  tuning,  and  development  on 
virtually  all  Boeing  commercial  jet  transports.  An  integral  part  of  its  use  has  been  a thorough 
verification  that  system  performance  on  the  computer  matches  that  obtained  during  air- 
plane flight  tests. 

Verification  of  simulator  performance  consists  of  comparing  results  with  results  from  similar 
tests  conducted  on  the  airplane.  In  this  study,  the  Boeing  airplane  performance  data  was 
readily  available  from  in-house  flight  testing.  Data  from  the  Rain  Tire  test  program  was 
provided  for  the  F-4.  Limited  data  on  the  C-141 A was  available  from  the  Combat  Traction  I 
program.  The  specific  test  conditions  used  for  correlation  are  listed  in  Table  6.  Conditions 
duplicating  airplane  weight,  CG  location,  brake  application  speed,  and  others  were  run  on 
the  simulator  and  compared  with  these  records.  Key  parameters  were  adjusted  until  the 
desired  level  of  correspondence  was  obtained.  The  following  items  were  considered  in  eval- 
uating the  correlation: 

• Stopping  distance 

• Skidding  pressure 

• Number  of  skids 

• Depth  of  skids 
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• Rate  in  and  out  of  skids 


• General  control 

In  addition,  strut  fore  and  aft  and  airplane  pitching  freqeuncies  and  damping  ratios  were 
checked  when  airplane  data  was  available. 

No  attempt  was  made  to  exactly  duplicate  airplane  stopping  distance  or  measured  parameter 
traces.  Instead,  emphasis  was  placed  on  producing  the  same  control  characteristics.  Stopping 
distances  were  matched  by  adjusting  the  tire-to-ground  friction  coefficient.  The  other  para- 
meters were  then  checked  and,  if  necessary,  adjusted  to  the  levels  measured  from  flight  test 
records.  Correct  skidding  pressure  levels  are  important  because  they  relate  brake  characteris- 
tics, runway  friction  levels,  and  tire  vertical  loading.  Similarly,  data  on  the  depth,  rates  in 
and  out,  and  the  number  of  skids  ensure  that  the  tire  simulation  is  behaving  as  it  should  and 
that  the  ground  force  simulation  is  properly  adjusted  for  the  particular  airplane. 

In  general,  traces  from  the  computer  tend  to  be  regular  in  nature  and  consistent.  Flight  test 
records,  though,  show  the  random  effects  produced  by  runway  discontinuities,  irregularities, 
and  second-order  effects.  Overall,  the  simulator  reproduces  the  predominant  factors  in  the 
stop. 

The  following  paragraphs  briefly  describe  the  correlation  of  each  airplane. 

I.  727 

Portions  of  test  181-16  condition  3.15.49.2.5  are  shown  in  Figure  11  along  with  an  equiv- 
alent condition  from  the  simulator.  The  control  is  characterized  by  multiple  skids  at  the 
start  followed  by  skids  at  essentially  regular  intervals.  The  skid  rate  is  about  once  per 
second.  Longer  intervals  occur  after  the  initial  cycling  (more  pronounced  on  some  wheels 
than  others).  Between  cycles,  a fairly  constant  pressure  increase  occurs  at  a rate  of  about 
250  psi/sec.  Skidding  pressures  are  about  1000  to  1 200  psi.  Both  the  flight  test  and  simula- 
tors show  these  characteristics.  The  airplane  braking  distance  was  1 163  ft,  and  the  distance 
obtained  from  the  simulator  was  1278  ft. 


2.  737 

The  737  simulator  performance  was  compared  against  Boeing  flight  test  29-35  condition 
1.20.05.001.1,  which  was  a performance  landing  at  a weight  of  100,500  lb.  The  center  of 
gravity  was  in  a forward  position  (5.3%  MAC),  and  brakes  were  applied  at  198  fps.  Portions 
of  the  wheel  velocity  and  brake  pressure  traces  from  this  test  are  shown  in  Figure  1 2 along 
with  corresponding  traces  from  the  computer. 

Analysis  of  the  records  showed  the  basic  similarity  between  the  airplane  and  the  simulator. 
Skidding  pressures  during  the  flight  test  ranged  from  2400  to  2800  psi  on  I he  left  side  of  the 
airplane,  while  those  on  the  right  side  ranged  between  1 800  and  2500  psi.  The  skidding 
pressures  on  the  simulator  ranged  from  2100  to  2400  psi.  At  high  speed,  the  skid  depths  on 
the  airplane  and  the  simulator  are  shallow  (1 5 to  20%  slip);  at  low  speed,  the  skids  become 
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much  deeper  (70-80%).  Rates  into  skids  range  from  400  to  600  fps-  and  rates  out  of  skids, 
from  600  to  1200  fps-  for  both  the  airplane  and  the  simulator.  Stopping  distances  were  also 
close,  with  1 387  ft  on  the  computer  compared  to  1 348  ft  for  the  airplane. 


The  747  correlation  runs  were  made  using  the  basic  747  antiskid  system.  This  system  uses  a 
50-tooth  wheel  speed  transducer;  the  sensitivity  tests  used  the  improved  system,  which  has  a 
200-tooth  transducer  and  incorporates  other  tuning  changes.  Extensive  use  of  the  simulator 
was  made  in  developing  the  improved  system  from  the  basic.  Detailed  correlation  of  the 
simulation  was  made  at  that  time. 

Foi  these  tests,  the  simulator  was  compared  against  two  flight  test  conditions.  Test  ! 1-8 
condition  1.20.052.003.5  was  a performance  landing  on  a dry  runway  with  an  airplane 
weight  of  486,100  lb.  Test  23-16  was  a stop  made  on  a wet  runway. 

For  Test  1 1-8  and  the  corresponding  simulator  run  (Figure  13j,  the  control  is  characterized 
by  multiple  skids  at  the  beginning.  Once  control  is  established,  skids  occur  at  essentially 
regular  intervals  about  I sec  apart.  At  high  speed,  the  skids  are  shallow,  becoming  deeper  at 
low  speed  even  to  the  point  of  temporary  lockups.  Likewise,  the  amount  of  pressure 
released  during  a skid  increases  until,  at  low  speed,  the  pressure  is  reduced  from  zero  during 
each  skid.  Skidding  pressures  are  around  1500  psi  at  high  speed,  increasing  to  around  2200 
psi  at  low  speed.  Stopping  distance  for  the  simulator  and  the  airplane  were  identical  at  2743 
ft. 

On  a wet  runway  (Figure  14),  the  wheels  begin  a series  of  deep  skids,  During  each  skid,  the 
brake  is  completely  released  and  a long  period  elapses  before  the  pressure  starts  to  increase. 
As  soon  as  the  pressure  starts  increasing,  the  wheel  goes  into  another  skid.  This  pattern  is 
repeated  during  most  of  the  stop. 


The  C-141A  brake  control  system  is  unique  among  those  tested  and  posed  unique  problems 
in  simulator  setup  and  correlation.  The  primary  concern  during  simulator  setup  was  t..~ 
wheel  speed  transducer.  Unlike  the  other  systems  tested,  the  C-141A  uses  an  AC  generator 
to  measure  wheel  velocity.  The  output  of  the  generator  is  an  AC  voltage  with  an  amplitude 
and  frequency  proportional  to  wheel  velocity.  The  output  is  rectified  and  filtered  in  the 
control  box  to  produce  a DC  voltage  proportional  to  wheel  velocity.  The  input  to  the  box 
from  the  computer  then  had  to  produce  the  same  amplitude  and  frequency  as  the  trans- 
ducer. (The  other  systems  worked  strictly  on  frequency  content  and  were  essentially  insensi- 
tive to  signal  amplitude.) 

To  ensure  correct  simulation  of  the  transducer,  an  aircraft  transducer  was  set  up  on  a test 
rig,  and  the  output  frequency  and  amplitude  were  determined  as  a function  of  velocity.  This 
output  was  then  produced  as  shown  in  Figure  1 5. 

Flight  test  records  for  the  C-141A  were  scarce. Those  available  were  from  the  Combat 
Traction  Flight  Test  program.  A dry  and  a wet  condition  are  shown  in  Figures  1 6 and  17. 
Because  of  the  scale  of  the  records,  little  quantitative  data  could  be  taken  from  them.  They 
do,  though,  provide  qualitative  data  on  the  general  control  of  the  system. 
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Figure  15.— C- 141  Wheel  Speed  Transducer  Simulation 
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Figure  16.-C-141  Dry  Runway 
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Figure  17.-C-141  Wet  Runway 
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Considering  the  dry  stop,  the  control  is  characterized  by  numerous  full  brake  dumps  (activa- 
tion of  the  Step  2 solenoid)  followed  by  a gradual  pressure  reapplication.  The  Step  1 
solenoid  (slower  partial  dump)  actuates  occasionally.  The  skidding  rate  is  about  one  skid  per 
second,  and  total  stopping  time  is  about  13  sec.  The  wheel  speed  traces  show  shallow  skids 
at  high  speed  and  deep  skids  at  low  speed.  The  actual  speed  departures  on  the  airplane  were 
probably  larger  than  the  traces  indicate,  because  the  low  transducer  frequency  requires  a 
large  amount  of  filtering. 

A dry  condition  run  on  the  simulator  is  also  shown  in  Figure  16.  As  in  the  flight  test  record, 
the  deep  dumps  predominate,  accompanied  by  an  occasional  activation  of  the  shallow 
dump.  Skidding  rate  is  about  one  skid  per  second.  The  wheel  speed  departures  are  much 
greater  than  those  from  the  flight  test,  as  expected.  The  trend  of  shallow  skids  at  high  speed 
and  complete  lockups  at  low  speed  is  duplicated.  The  stopping  distance  for  this  condition 
was  1841  ft,  compared  to  a handbook  value  of  1900  ft. 

The  wet  runway  landing  (Figure  17)  shows  a continuous  sequence  of  deep  skids  at  a rate  of 
about  three  skids  per  second.  The  control  is  characterized  by  a reapplication  pressure  higher 
than  the  runway  can  sustain.  Hence,  the  wheel  goes  into  a skid  as  soon  as  pressure  is 
reapplied.  Similarly,  the  simulator  records  show  the  same  dump-fill  pattern  at  about  the 
same  rate. 


5.  F-4 

Control  of  the  F-4  on  a dry  runway  (Figure  18)  is  characterized  by  a large  number  of 
shallow  -kids.  Most  of  the  skids  result  in  a short,  partial  dump  of  brake  pressure.  As  speed 
decreases,  the  skids  generally  become  deeper  and  fewer.  Another  prominent  characteristic  of 
this  airplane  is  the  gradual  increase  in  skidding  pressure  as  the  speed  decreases.  The  aircraft 
braking  distance  for  this  control  was  3400  ft,  while  the  simulator  distance  was  3278  ft. 

On  a wet  runway,  Run  23,  the  skids  are  much  more  definite,  being  deeper  and  of  longer 
duration  (Figure  19).  Only  one  wheel  on  the  airplane  goes  into  a skid.  This  is  caused  by  the 
paired  wheel  control  used  on  the  airplane.'  The  brake  pressure  is  operating  in  a dump-fill 
mode  common  to  many  systems  on  low  mu  surfaces.  Equivalent  conditions  from  the  simula- 
tor show  the  same  characteristics.  The  aircraft  braking  distance  was  6221  ft,  and  the 
simulator  distance  was  5890  ft. 

6.  OPERATIONAL  ANOMALIES 

The  basic  correlation  described  in  preceding  paragraphs  has  been  limited  to  the  normal 
operation  of  brake  control  systems.  Under  certain  conditions,  a system  does  not  operate 
normally.  These  anomalies  exist  only  under  extreme  operating  conditions. 

Two  anomalies  will  be  discussed  here  as  a further  indication  of  the  ability  of  the  simulator 
to  duplicate  the  behavior  of  the  airplane. 
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The  first  anomaly  involves  the  operation  of  the  locked-wheel  protection  device  in  the  727 
Mark  II  antiskid  system.  Locked-wheel  protection  normally  prevents  complete  wheel  lock- 
ups. But  if  the  pilot,  on  a runway  with  very  low  friction,  meters  full  pressure  to  the  brakes 
and  the  two  paired  wheels  spin  down  together,  the  locked-wheel  protection  for  those  wheels 
can  disarm  and  the  wheels  can  lock.  This  was  demonstrated  during  USAF/NASA/FAA 
testing  at  Edwards  AFB.  Figure  20  shows  wheel  velocities  and  brake  pressures  for  the 
airplane  and  similar  conditions  on  the  simulator. 

The  second  anomaly  relates  to  a 737  landing  at  Rosewell,  New  Mexico,  on  a wet  runway. 
The  brakes  were  applied  before  full  wheel  spinup,  causing  the  system  to  operate  for  an 
extended  period  of  time  in  the  high-slip  region.  Wheel  velocity  and  brake  pressure  traces  for 
this  condition  and  a computer  run  are  shown  in  Figure  21 . 

While  these  two  cases  show  that  the  simulator  will  reproduce  abnormal  conditions,  this 
study  was  primarily  concerned  with  normal  system  operation.  Abnormal  conditions  must  be 
considered  case  by  case  because  factors  that  may  otherwise  have  only  minor  influence  can 
become  major. 


SECTION  IX 

DETERMINATION  OF  SIGNIFICANT  PARAMETERS 


1.  PARAMETER  RATINGS 

The  parameter  rating  index  (PR1)  was  used  to  rank  the  parameter  changes  according  to  their 
effect  on  braking  distance.  The  value  of  the  PRI  is  the  average  percentage  deviation  from  the 
baseline  braking  distance.  Thus,  the  larger  the  PRI,  the  greater  impact  the  parameter  has  on 
airplane  braking  distance.  Based  on  the  PRI,  the  parameter  changes  have  been  arranged  in 
numerical  order,  and  the  results  are  listed  in  Table  7.  The  table  presents  the  average  PRi  for 
the  five  aircraft  tested.  The  PRI  for  each  airplane  taken  individually  is  listed  in  ASD- 
TR-74-41,  Volume  II,  Section  XII.  Although  the  tire-ground  friction  coefficient  (mu)  is  a 
predominant  influence  on  airplane  braking  performance  (see  Figure  9),  it  is  not  included  in 
Table  7 because  the  PRI  rating  methodology  was  used  to  determine  variable  significance  for 
a range  of  mu  values. 
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Table  7.- Parameter  Rating 


Description 
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2.  SIGNIFICANT  PARAMETERS 

The  second  step  in  the  rating  of  parameters  was  to  determine  which  parameter  changes  have 
a significant  effect  on  stopping  distance.  It  was  decided  to  consider  all  parameters  having  a 
PRI  greater  than  2.0  as  being  significant.  A value  of  2.0  represents  a 2%  change  in  the 
baseline  braking  distance.  The  repeatabilit-  e analog-hardware  simulation  itself  results 
in  \%  variations. 

Based  on  the  above  criterion,  Table  7 has  bee,,  cduici1  , cak  available  mu  included.  In 
addition,  the  parameter  changes  have  been  summari  «d  ' ••  nbining  related  tests  under  a 

general  heading.  The  resulting  list  of  parameters  havi  lificant  effect  on  braking 

distance  is  given  in  Table  8. 

The  directional  control  impact  of  some  of  these  variables  was  also  assessed,  but  the  data  is 
not  shown  in  the  form  of  bar  charts.  The  indices  of  directional  performance  measured 
during  this  program,  although  significant,  are  not  the  only  contributing  factors.  Undue 
emphasis  on  these  could  lead  to  misunderstanding  rather  than  clarification. 

Generally,  inefficient  skid  control  systems  display  high  skidding  values,  thus  reducing 
cornering  capability;  a decrease  in  cornering  coefficient  significantly  degrades  directional 
control.  True  assessment  of  this  aspect  will  require  an  in-depth  study  using  steering  and 
ground  handling  simulators.  The  reader  is  therefore  cautioned  regarding  the  use  of  cornering 
index  as  noted  during  the  tests.  See  ASD-TR-74-41,  Volume  II,  Section  XII,  for  numerical 
values  of  the  cornering  index  for  all  aircraft  tested. 


Table  8. -Reduced  Significant  Parameters 


C 

Punk  available  ground  friction 

2, 

Spoiler  or  drag  device  effectiveness 

3. 

Haad  or  tailwind 

4. 

Brake  application  speed 

5. 

Brake  pressure  application  rate 

6, 

Landing  weight 

7. 

Metered  pressure  effectiveness 

8. 

Drag  device  deployment  timing 

9, 

Mu -slip  curve  shapes: 
Tire  In ilation  pressure 
Tire  hooting 
Flat  /i-r/penk 
Peak  fx  percentage 

10. 

Ambient  temporature 

11. 

Engine  Idle  thrust 

12, 

Contar  of  gravity  location 

13, 

Brake  pressure  application  timing 

...  J 
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SECTION  X 

SELECTION  OF  PERTINENT  PARAMETERS 


Figure  22  is  a flow  chart  where  each  block  represents  a major  step  of  analysis  in  the 
formulation  of  the  prediction  equation  and  the  resulting  methodology. 


Table  9 lists  significant  parameters  after  excluding  the  parameters  that  depend  on  pilot 
technique  and  are  outside  the  scope  of  the  present  work.  The  following  paragraphs  present 
the  reasoning  for  further  refinement  to  the  list. 


1 

1 


Table  9.— Significant  Parameters  * 

1.  Peak  available  mu 

2.  Drag  device  offectlvenesi 

3.  Head  or  tall  wind 

4.  Brake  application  speed 

6.  Landing  weight 

6,  Air  density 

7,  Engine  idle  thrust 

8,  Center  ot -gravity  location 


* Excludes  parameters  dependent  on 
pilot  technique. 


1.  BRAKE  APPLICATION  VELOCITY 


If  an  airplane  is  loaded  heavily,  its  landing  speed  increases.  Because  the  maximum  lift 
coefficient  and  wing  urea  remain  the  sume,  ‘he  landing  speed  Is  a direct  function  of  gross 
weight.  Accordingly,  for  u lower  landing  weight,  the  landing  speed  decreases.  So  the  landing 
weight  of  an  airplane  cannot  he  varied  independent  of  the  landing  velocity.  Because  pilot 
techniques  vary,  however,  the  velocity  can  vary  independent  of  weight. 

One  of  the  requirements  of  forming  dimensionless  groups  is  that  each  term  or  group  be 
independent.  Therefore,  only  velocity  (but  not  weight  und  velocity)  wus  chosen  as  a signifi- 
cant independent  variable. 


2.  AERODYNAMIC  LIFT  AND  DRAG 

The  same  reasoning  applies  to  the  variation  of  lift  and  drag  coefficients.  One  cannot  be 
varied  without  varying  the  other,  so  the  C'lK'd  ratio  was  chosen  us  un  independent  variable. 
The  only  exception  If  the  F-4  where  no  spoilers  und  only  a drag  chute  Is  employed  for 
uerodynamlc  braking.  The  C'l  thus  remains  constant.  However,  when  C|j  was  considered  is 
un  Independent  variable  tirterm)  instead  of  ('l/C|)  lor  the  F-4,  the  predlcltloii  equation  for 
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the  model  did  not  change  in  its  accuracy  or  complexity.  Thus,  the  Cj_/C  q it  term  was  used 
for  F-4,  as  for  the  other  airplanes  studied. 

3.  WIND  VELOCITY 

Head  and  tail  winds  affect  the  air  distance  and  transition  distance  as  well  as  the  braking 
distance.  The  analog  simulation  and  dimensional  analysis  represent  only  the  braking  segment 
of  the  landing  process.  The  question  therefore  was  whether  the  wind  velocity  should  be 
considered  as  an  independent  variable  for  the  dimensional  analysis. 

The  effect  of  brake  application  speed  and  that  of  the  head  and  tail  winds  on  the  stopping 
distance  was  plotted  and  compared  (Figures  23  and  24).  In  Figure  23,  for  example,  curves 
have  been  drawn  through  737  and  C-141  data  points;  similar  curves  could  be  drawn  for  the 
other  airplane  data  points.  Since  a single  curve  can  describe  the  variation  of  both  wind 
velocity  and  brake  application  speed,  wind  velocity  need  not  be  considered  as  an  independ- 
ent variable. 


4.  CENTER-OF-GRAVITY 


The  CG  location  effect  considered  in  the  analog  simulation  was  of  the  form: 


\]f  B 


m*la 

LA  + LB+g-HB 


where  LA,  -LB,  and  HB  are  geometric  distances  (see  Table  2),  and  p is  the  coefficient  of 
available  friction.  Because  the  CG  variation  of  any  correlation  ground  vehicle  would  be 
minimal,  and  because  geometric  similarity  with  an  airplane  would  be  almost  impossible  to 
achieve,  it  was  decided  to  consider  only  the  coefficient  of  available  friction  as  the 
independent  variable. 


5.  OTHER  PARAMETERS 

Tire  remaining  parameters  (p  and  F0)  are  independent  variables  and  require  no  discussion. 
From  the  preceding  paragraphs,  it  follows  that  the  pertinent  variables  are: 


Brak'ng  disi  nice 

(s) 

Available  mu 

<p) 

CL/CU  ratio 

(CL/CD> 

Brake  application  speed 

(v) 

Air  density 

<p) 

Engine  idle  thrust 

(F'c) 

Velocity,  V (knots) 
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SECTION  XI 

DEVELOPMENT  OF  PREDICTION  MODEL 


1.  COMPONENT  EQUATIONS 

The  first  step  in  forming  a prediction  model  is  to  identify  the  pertinent  and  independent 
variables.  This  step  is  by  far  the  most  important  because  the  validity  of  the  results  depends 
on  the  correctness  with  which  the  pertinent  factors  are  selected.  For  this  study,  as  explained 
in  Section  X,  this  required  ht  combining  of  some  of  the  interdependent  variables  listed  in 
Table  9.  The  list  of  resultan,  independent  variables  is  shown  in  Table  10. 

The  second  step  is  to  express  the  secondary  quantity  (dependent  variable)  as  a function  of 
the  primary  quantities  (independent  variables),  so  that: 

s = F (g,  v,  p,  Fe,  ft,  C^/Cp)  ( i ) 


where: 


s = braking  stop  distance 
g « acceleration  cuused  by  gravity 

The  dimensional  matrix  that  can  be  formed  for  the  fundamental  units  (mass,  length,  and 
time)  of  the  seven  parameters  in  Eq  1 is  of  rank  3,  so  that,  according  to  Buckingham’s  ir 
theorem,  these  would  yield  four  independent  it  terms.  By  inspection  and  analysis,  they  can 
be  written  (sg/v2),  (p),  (Cl/C’d>,  and  (pv(’/Feg2).  Thus: 

(sg/v2)  = F (p,  CL/CD,  pv^F,,  g2) 
or:  (it  | ) = F ^ir2,  Itq,  ir^ 

where: 

IT  | = Sg/V2 

ir2  = g 

ff3  = CL/CD 
*4  “ pv6/ Fe  g2 

Appendix  C shows  the  detailed  unalysis  of  arriving  at  Eq  2 and  Eq  2u.  The  application  of 
dimensional  anulysis,  including  tile  pi  theorem,  leads  to  a type  of  equation  involving  an 
unknown  function,  of  which  Eq  2 is  un  example.  Before  a prediction  equation  can  be 
formulated,  the  nature  of  the  function  must  be  determined.  This  cunnot  be  accomplished  by 
dimensional  anulysis,  but  it  can  be  done  from  anulysis  of  laboratory  obscivutions. 


(2) 

(2a) 
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Table  lO.-Pertinent  Independent  Variables 


Variable 

Notation 

Available  mu 

M 

Brake  application  speed 

V 

Drag  device  effectiveness 

CL/CD 

Engine  thrust .... 

Fe 

Air  density 

P 
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The  best  procedure  for  evaluating  a function  is  to  arrange  the  observations  so  that  all  but 
one  of  the  pi  ter.ns  containing  the  independent  variables  in  the  function  remain  constant. 
Then  the  remaining  independently  variable  pi  term  is  varied  to  establish  a relationship 
between  it  and  the  dependent  variable  (jt|  term).  This  procedure  is  repeated  for  each  of  the 
pi  terms  in  the  function;  the  resulting  relationships  between  ir|  and  the  other  individual  pi 
terms  are  called  component  equations.  Statistical  curve  fitting  computer  programs  were  used 
to  generate  the  component  equations  (see  Appendix  D)  A summary  of  the  equations  is 
listed  in  Table  1 1. 


2.  GENERALIZED  FUNCTIONS 

When  the  component  equations  have  been  determined,  they  are  combined  in  a certain 
manner  to  give  a general  relationship,  It  is  possible  for  some  of  the  component  equations  to 
be  combined  by  multiplication,  while  others  require  addition  in  the  formation  of  the 
resultant  prediction  equation.  In  general,  these  two  methods  are  adequate  for  the  majority 
of  engineering  problems.  For  the  stopping  distance  problem,  the  analysis  showed  that  the 
prediction  equation  should  be  formed  by  multiplication.  The  necessary  and  sufficient  condi- 
tions to  be  met  for  the  function  to  be  a product  were  developed  and  translated  into  tests  of 
validity.  All  aspects  of  the  development  of  prediction  equations  discussed  in  this  paragraph 
are  detailed  in  Appendix  E.  The  major  equations  of  interest  are  repeated  in  succeeding 
paragraphs. 

When  the  component  equations  (see  Table  1 1)  are  combined  by  multiplication,  the  predic- 
tion equation  is  of  the  form; 

»|,  = (CI  I ' 3.4  (,r  I ' 2.4  <",)  2,3  (40) 

where  tire  bar  denotes  a constant  (held)  value, 


Tlte  analysis  shows  that  the  value  of  the  constant  term  C is  of  tire  form 


[F  (f2,  rrp  rr4)]  2 

Thus  tile  prediction  equation  is  of  the  form; 


F (ir2,  *3,  ir4) 


F (1r2.rr3.Tf4)  F (f2,  ir3,ir4)  F (i2,  *3,  rr4) 
[F  (ff2.  rr3.  rr4)]  2 


(41) 


(42) 
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Airplane 
model  p * 


737  0.6 


747  0.6 


Table  1 1. -Summary  of  Component  Equations 


Or, I -0.7048  0r2)  -0'8196 

(ff,l  - 1,06869  (ir3l  10  312B-0.33’B  % SP]  ** 

Or,)  = 2.314  (ir4l [-0.06836] 

Or,) 1 0.7048  0r2)(0-a196! 

Iff,)-  1.4882  0r3)  (0.3090.0.3070%  SP| 

Iff,)-  4.2724  (ff4)'0'00409 

Iff,)  » 0.6648  (ff2l  -1'125 

Iff,) - 3.7262  lff3l  (0.21686.0.13473  % SPl 

Iff,) - 27.287  Ola)  amB 

Iff, )- 0.771 6 0r2)  °'7647 

(ff,l-  1.1693  <ff3)  I0.168B3-0.0181 13  %SP| 

Iff,  I - 1.4078  0r4)  -0.01951 


Or,) -0.7718  Or,)  °'7947 

Iff,)  * 1,6637  (ffg)  (0  1 8666-0. 0B8B3  % SP | 

Iff,)-  2.1768  (ff4l  -0-0327 

Iff, I « 0.7716  0r2)  0 2647 

Iff,)-  2.868/ lff3)  (0  18397.0.06279%  SP| 

01,1-  4.8621  lff4)  0,08739 


Or,  I - 0,838  0r2)  u'81s 

Iff,)-  1.2214  (ff3l  10, 14803-0, 1 1286%  SP| 

Or,)  - 1.6182  Or4l  -0-01676 


Value  of  p used  In  the  data  set. 

% SP  Is  the  percentage  of  the  spoiler  configuration, 


Table  1 1. -Summary  of  Component  Equations  ( Continued ) 


Equation 

Eq 

No. 

(7T,)  - 0.838  0r2r0-81 6 

(18) 

(*,>- 1.5641  ,#3)  10.1 8957.0.091 83  % SP] 

(21) 

(IT,)  = 2,3966  <jr4)‘0  03004 

(22) 

(jr,)  - 0.838  (jr2>  0 816 

(18) 

(7T,)  - 2.6185  (,T3)  t0.2365-0.0818  % OP) 

(23) 

01,1  -5.4099  Or4)°’06B88 

(24) 

Or,)  -0.876  Orj)'1  0268 

(26) 

(7T,)  - 1.402  (W3>  [0.2134^.13084  % SP] 

(28) 

01,1  - 2.8836  0r4)  0 06,29 

(27) 

Or,)  = 0.876  0r2r1  0288 

Or,)-  2.1846  (ff3)I0.23065.0.14819  96SP) 

IKiH 

(291 

Or, 1 =6.0891  (rr4)  ■°8694 

Or,)  = 0.876  Or/1  0288 

(26) 

(ir, ) = 4,51737  (rr3)  (0.21 61 4-0.1 8567  % SPJ 

(30) 

(ir,)  = 21.6248  (jt4)  0-14344 

(31) 

Or,)  =0.84576  <7r2)°  9239 

(32) 

Or,)  = 1.4699  Or3)+0'3036 

(33) 

Or,)  =2.5323  (7T4)004408 

(34) 

Or,)  = 0.84575  0r2)'°’9239 

(32) 

<7T,)  - 2.1616  (jr3!  l0.52B98-+O.Oa287  % SP] 

":5i 

(w,)  - 5.8872  (ir4)  00789B 

(36! 

The  equations  constituting  a test  for  the  validity  of  Eq  42  are  shown  to  be  (see  Appendix  E): 


F (if2,  tr3,  tr4)  F (vf2,  rr3,  rr4)  F (ff2,  ir3,lr4)  F 4>  (43) 

[F  (jr2,?3,W4)]2  “ [F  (?2,tr3.W4)]2 

or: 

F (ir2,  if3,  ?4)  F (ff2,  ff3,  tr4)  F 0r2,#3,tr4)  F (f2,f3,  ir4) 
= (43a) 

[F  *3*  ?4>]  2 [F  <*2'  "3.  »4>]  2 

The  values  5r2  and  tr3  are  values  of  ir2  and  tr3  held  constant  at  some  value  other  than  ?->  and 
W3.  Thus  from  the  observed  data: 

ir2  » 0.6  the  primary  set  of  data,  for  example 

#2  » 0.4) 

| Q V supplementary  sets  of  data 

If  the  supplementary  sets  of  data  satisfy  either  43  or  43a,  the  general  equation  can  be 
formed  by  multiplying  the  component  equations  together  and  dividing  by  the  constant,  as 
indicated  in  Eq  42. 

This  test  was  applied  to  all  available  data  (component  equations);  the  results  are  shown  in 
table  12,  clearly  indicating  the  validity  of  the  approach.  Table  E-l  contains  the  details  of 
this  calculation. 

Another  test  of  validity  was  to  calculate  the  value  of  the  constant  term  C of  Eq  40.  The  test 
requires  that  any  of  the  three  component  equations  should  yield  an  identical  value  for  C. 
This  test  was  also  applied  to  all  the  test  data;  the  results  are  showr  - Table  13.  Again,  the 
accuracy  achieved  is  satisfactory.  Table  E-2  contains  details  of  this  calculation. 

The  two  validity  tests  were  successful,  thus  permitting  the  writing  of  the  prediction 
equations.  A summary  of  all  prediction  equations  is  listed  in  Table  14.  Equation  44  is  a 
combination  of  Eq  3,  4,  5 and  corresponding  C.  Equation  45  is  a c>  bination  of  Eq  3,  6,  7 
and  corresponding  C,  and  so  on. 
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Table  14.—  Summary  of  Prediction  Equations 


Airplane 

model 

n' 

Equation 

Ml 

727 

0.6 

hr,)  = 1.5001  <7r2)-°-8196  ,73)i°'3128'0-3375  % SP]  ,^0.06836 

(44) 

m 

Or,)  = 2.0197  (7j)0fi'98  („3,t0.3090.0.30107%SP)(l4)  0.0941 

(45) 

(7,)  - 4.558  (72)'-126  ,ff3,[O.2159.0.1347  %SP)  ,^-0.1768 

(46) 

737 

0.6 

(77,1-0.9632  (tt2)  0 7647  ^ 10.1567-0.081 13  % SP  | ,^,-0.01 951 

(47) 

(*,)  . 1.096  !ir2)-°'7647  (ir3)  [0.16655.0.08663  % SP] , ^,  0.03222 

(481 

rii 

(77,1  . 1.4409  (7T2)  0 7647  (ff3!t0.16397.0.06279%SPl  ^fl.06739 

(49) 

747 

0.6 

(7,1  - 0.9539  (772)  0 815  (73) E0. 14803-0. 11 286  %SP1  ^ -0.01676 

(50) 

0.4 

(*,)  - 0.9974  <72>°-815  (ir3l  [0.18957-0.09193  % SP1  ,^,0. 03004 

(51) 

0.2 

(77,1  = 1.1897  (72)  0'816  (tt31  {0.2365-0  0818  % SP1  ,^,  0.05688 

(52) 

C-141 

0.6 

(77, 1 * 1.6814  (ff3,t0.21341-0.13084%SP)  ,^,0.06129 

(53) 

0.4 

(7,1  - 2.20116  (772)  1 0268  10.23065-0  1 48 1 9 % SP)  ,^,-0.08694 

(54) 

0.2 

(77,1  - 4.1160  (J72l1  0268  (773)  [0.2 1 01 4-0. 1 6567  % SP]  ,^,-0.14344 

155) 

F-4 

0.6 

(17,1  = 1.6/90  (7721  0 9239  (773)  0 3636  (774)  0 04408 

156) 

0.4 

(7,1  - 2.80256  (772)  O-9239  (773) [0.52698+0.08287  %SP1  ^i'007895 

(57) 

0.2 

(7,1  = 31.52  (72)'1  0694  |73l  0'69738  (J74r°  26017 

(58) 

* Value  used  to  derive  the  equation. 
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SECTION  XII 

MODEL-TO-SIMULATOR  CORRELATION 


The  prediction  equations  were  next  used  to  correlate  back  with  the  stopping  distance  data 
collected  in  the  Task  I simulation.  A summary  of  errors  in  correlation  is  listed  in  Table  1 5. 

The  limitations  (range  of  validity)  of  the  prediction  equations  are: 

• Equations  47  through  55  are  applicable  for  p values  of  0. 1 to  0.6. 

• Equations  44,  45,  56,  and  57  are  applicable  for  p values  of  0.3  to  0.6 

• Equations  4 > and  58  are  applicable  for  p values  of  0. 1 to  0.2  only. 

For  a given  airplane  model,  the  three  prediction  equations  are  interchangeable,  alternate 
solutions  if  their  range  of  applicability  and  validity  is  common.  Thus,  Eq  45  and  58  are 
unique  solutions  and  not  interchangeable  with  their  counterpart  Eq  44  and  4S,  or  56  and 
57.  Some  airplane  systems  need  only  one  prediction  equation  to  define  the  entire  range  of  p 
values  tested  on  the  simulator;  others  needed  more  than  one  equation.  The  reason  for  this 
can  be  comprehended  by  studying  braking  distance  efficiency  curves  for  the  various  systems 
as  shown  by  Figure  25. 

Braking  distance  efficiency,  r)s,  is  defined  as  the  ratio  or  the  perfect  braking  distance  to  the 
braked  airplane  distance  resulting  from  the  simulation. 

rjs  = Xp/Xa  x 100% 


rjs  = braking  distance  efficiency 
Xp  = perfect  braking  distance 
Xa  = airplane  braking  distance 


The  perfect  braking  distance  is  the  distance  required  to  stop  the  airplane  if  it  is  braked  for 
the  entire  stop  with  maximum  available  braking  force.  Braking  distance  efficiency  indicates 
the  degree  to  which  the  system  meets  its  primary  requirement  of  stopping  the  aircraft. 

The  skid  control  systems  for  the  five  subject  airplanes  encompass  three  generations  of 
technology  namely  the  old,  intermediate,  and  an  advanced  type.  A study  of  Figure  25  shows 
that  the  curves  for  the  advanced  technoiogy  system  is  nearly  linear,  sharp  changes  in  slope 
appear  in  the  curves  for  the  old  as  well  as  the  intermediate  technology  systems  at  p values  of 
0.3.  Piecemeal  linearization  is  always  required  when  writing  mathematical  relationships  for 
curves  of  the  type  shown  for  the  old  and  intermediate  technology  systems.  That  is  why 
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data  at 

i%) 

727 

(44) 

0.8 

•2.9  to  +4.8 

0.4 

•5.0  to  +3.5 

(45) 

0.4 

-4.8  to  +1.9 

0.6 

•4.9  to  +5.0 

(46! 

0.2 

•4.9  to  +3.7 

737 

(47) 

0.6 

-1.7  to +1.8 

0.4 

-1.2  to  +3.9 

0.2 

-3.6  to  +3.1 

(48) 

0.6 

-3.4  to  +1.3 

0.4 

• 1.0  to +1.8 

0.2 

-3.1  to  +4.6 

(48! 

0.6 

■4.9  to +1.8 

0.4 

•2.9  to  +2.6 

0.2 

■2.0  to  +2.0 

(56) 

0.6 

•2.9  to  +0.7 

0.4 

-5.6  to  +5.0 

(57) 

0.4 

•5.4  to  +5.4 

0.6 

-6.7  to  +7.3 

(68) 

0.2 

-1.7  to  +3.3 
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Peak  available  Mu 


Figure  25.  -Mu-E fficiency  Curves 


more  than  one  prediction  equation  was  necessary  to  treat  these  two  cases.  Performance 
deficiency  of  the  old  and  intermediate  systems  at  low  u values  has  led  to  the  development  of 
advanced  systems. 

The  correlation  data  error  summary  (Table  1 5)  indicates  that,  for  almost  all  conditions,  a 
prediction  accuracy  of  ±5%  can  be  achieved. 

Even  though  in  the  correlation  process,  comparison  was  made  between  predicted  and  actual 
ttj  values,  that  is: 

*)  •) 

(sg/v")pred  vs  (sg/v-)actual 

it  is  tantamount  to  comparing  the  braking  stop  distances  because  both  terms  use  identical  g 
and  v values,  and  the  distance  term  s has  no  exponent.  For  example,  from  727  correlation 
calculations: 


(ff  | )pred  = 1°78;  (it, )acluai  = 1 .082;  error  percentage  = -0.4 
(s)pred  = 1273  ft;  (s)actua|  = 1 278  ft;  error  percentage  = -0.4 
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SECTION  XIII 

MODEL-TO-SIMULATOR  QUALITATIVE  COMPARISON 


I.  PREDICTED  MU  FROM  FLIGHT  TEST  DATA 

The  credibility  ol'  the  simulator  was  established  by  comparing  the  simulator  and  airplane 
flight  test  data  and  showing  that  similar  trends  were  obtained  under  identical  conditions  (see 
Section  VIII).  The  credibility  of  the  prediction  model  has  been  established  by  obtaining  a 
±5%  correlation  accuracy  in  predicting  simulator  stopping  distances.  The  next  logical  step  is 
to  determine  if  the  airplane  flight  test  data  could  be  correlated  to  the  prediction  model.  The 
results  of  this  exercise  are  shown  in  Table  16.  From  the  type  of  information  available  on  the 
flight  test  data  (Table  16)  the  only  parameter  that  could  be  calculated  by  the  prediction 
equation  was  the  friction  coefficient.  The  predicted  values,  for  both  dry  and  wet  runway 
conditions,  appear  reasonable. 


2.  WET-RUNWAY  ANALYSIS 

During  Task  1 simulation  testing,  a wet  runway  was  simulated  so  that  the  available  ground 
mu  was  programmed  to  vary  with  speed.  (See  Figure  26).  The  mu  values  (end  points)  used 
were  0.05  at  brake  application  speed  and  0.5  at  the  end  of  the  stop. 

The  average  value  of  peak  available  mu  for  the  braking  system  was  unknown,  so  it  was 
decided  to  use  the  component  equations  formed  earlier  (the  tt|  versus  rrj  relationships)  to 
calculate  peak  available  mu.  Based  on  calculations  for  wet  runways,  prediction  equations 
were  generated  for  wet  runway  cases.  With  these  prediction  equations,  a correlation  predic- 
tion accuracy  analysis  was  conducted  as  before  and  satisfactory  results  were  obtained.  The 
details  of  this  analysis  are  reported  in  ASD-TR-74-41, Volume  II,  Section  XVI.  The  results 
gave  additional  confidence  to  the  selected  methodology  for  forming  prediction  equations. 


3.  DATA  VALIDITY 

Yet  another  test  of  the  prediction  equations  came  about  accidently.  During  the  correlation 
process,  in  spite  of  excellent  correlation  an  unusually  high  error  would  result  for  one  or  two 
conditions,  even  though  there  was  close  correlation  between  predicted  and  actual  distances, 
This  happened  with  three  airplanes.  When  the  raw  data  from  Task  I was  rechecked,  book- 
keeping errors  were  discovered  in  the  con  u lions  that  produced  the  anomalies.  Thus  the 
prediction  equation  was  aole  to  sort  out  the  unusual  cases  that  failed  to  fit  the  pattern, 
pointing  out  bad  data. 
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SECTION  XIV 

DETERMINATION  OF  VEHICLE  CRITERIA 


The  formulation  of  the  prediction  equation  has  been  accomplished  with  the  use  of  dimen- 
sional analysis.  In  arriving  at  the  equation,  a complex  dynamic  process  was  duplicateu  with  a 
hardware-analog  brake  control  simulation.  Comparison  between  actual  airplane  stops  and 
the  simulation  indicated  similar  braking  performance  and  established  a high  degree  of  con- 
fidence in  this  approach.  The  effect  of  parameter  variations  cn  stopping  distance  was  then 
examined  on  the  brake  control  simulator  and  compared  to  baseline  data.  Parameters  that 
changed  the  stopping  distance  by  more  than  2%  were  considered  significant  and  included  in 
the  dimensional  analysis.  The  resulting  prediction  equation  appears  with  the  general  format: 

sg/v2  = C(M)a(CL/CD)^(pv6/^g2)8 

The  equation  has  four  different  nondimensional  terms.  The  dependent  term  represents  a 
distance,  s,  which  is  normalized  with  the  gravitational  constant,  g,  and  the  aircraft  velocity, 
v,  relative  to  the  ground  at  brake  application.  Head  or  tailwinds  can  affect  the  airplane 
velocity  at  touchdown  and  must  be  taken  into  consideration,  in  addition,  pilot  technique 
plays  a major  part  in  the  determination  of  touchdown  velocity.  The  dimensionless  coeffi- 
cient, p,  is  the  peak  friction  efficient  generated  at  the  tire-runway  interface  under 
dynamic  conditions  of  the  brak.ng  process.  For  wet  runway  conditions,  this  "mu"  can  be 
obtained  using  the  technique  shown  in  ASD-TR-74-41 , Volume  11,  Section  XVI.  The  aero- 
dynamic term,  Cl/C[),  is  a ratio  of  lift  coefficient,  Cl,  and  drag  coefficient,  Cp,  of  the 
airplane  on  the  ground  in  braking  configuration.  Variations  in  these  coefficients  result  from 
the  use  of  different  flap  settings,  spoilers,  anc  .'rag  devices.  The  last  term,  pv*/ Feg2,  is  a 
combination  of  runway,  altitude,  and  temperature  information  as  expressed  in  air  density,  p, 
and  total  engine  idle  thrust  during  braking.  The  exponents  o,  (3  and  6 and  the  coefficient  C 
result  from  the  sensitivity  study  and  the  dimensitnal  analysis.  For  each  airplane,  a unique 
coefficient  and  set  of  exponents  exist. 

The  equation  permits  the  calculation  of  the  airplane  stopping  distance,  assuming  proper 
information  of  airplane  and  weather  parameters  and  an  accurate  and  meaningful  measure- 
ment or  prediction  of  the  tire-runway  friction  coefficient.  In  turn,  when  the  length  of  a 
specific  runway  is  known,  the  safety  margin  for  a landing  and  the  last  point  for  safe  brake 
application  can  be  determined.  The  most  elusive  value  in  the  equation  is  an  adequate  value 
for  the  tire-runway  friction  coefficient.  This  is  the  peak  value  developed  under  the  dynamic 
conditions  of  braking  and  is  represented  in  the  hardware-analog  simulation  with  a friction 
slip  curve.  The  dimensionless  appearance  of  the  prediction  equation  of  similtude  and  the 
friction  coefficient  must  be  derived  and  measured  for  a similar  process.  Hence  the  following 
conditions  must  be  observed. 

Deceleration  Process:  Tile  tire  is  deformed  by  shear  forces  in  the  braking  process,  and  a 
shift  in  the  footprint  occurs. 

Rotating  Wheel:  During  the  braking  process,  the  wheel  is  controlled  with  an  antiskid  device 
that  protects  it  from  lockup  and  continuously  hunts  for  the  peak  friction  value. 
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Energy  Absorbed  in  Brake:  This  implies  that  the  major  portion  of  the  energy  absorption  is 
taking  place  in  the  brake.  While  braking  can  be  realized  with  a locked  wheel,  the  true 
dynamic  process  must  be  similar  to  that  taking  place  on  the  braked  airplane  wheel. 

Distance  Generated:  A measurable  ground  stopping  distance  must  be  available  as  a data 
point.  Hence,  braking  should  be  continuous  from  the  time  of  brake  application  until  a full 
stop  is  reached. 

Aircraft  Tire:  While  similar  in  fun  :tion,  aircraft  and  automobile  tires  are  different  in 
construction.  Notably,  the  inflation  pressure  of  aircraft  tires  is  much  higher  than  in  auto- 
motive tires,  and  the  theoretical  hydroplaning  speed  is  directly  affected.  The  footprint  area 
of  an  automobile  tire  is  rectangular  with  a constant  width,  while  that  of  the  airplane  tire  is 
elliptical. 

Design  criteria  for  a ground  vehicle  used  in  the  prediction  process  must  consider  the  above 
factors.  Dynamic  siiniltude  could  be  obtained  with  the  use  of  a model  airplane  with  a brake 
system  similar  to  those  of  existing  systems.  This  would  result  in  a complex  and  costly 
approach;  in  many  cases,  it  would  not  be  practical  to  fly  the  model.  It  appears  that  the  most 
desirable  approach  would  involve  a vehicle  that  demonstrates  dynamic  similarity  in  the 
braking  process.  This  could  be  a ground  vehicle  equipped  with  airplane  tires  (Type  VII)  and 
a suitable  antiskid  system.  The  vehicle  could  be  used  in  predicting  a meaningful  tire-runway 
friction  coefficient  as  shown  in  Figure  27.  The  approach  outlined  here  is  similar  to  that  used 
in  arriving  at  the  airplane  prediction  equation.  Simulation  of  the  vehicle  and  its  brake 
control  system  would  permit  a sensitivity  study  to  define  a stopping  distance  prediction 
equation  based  on  the  most  important  parameters.  Most  likely,  the  equation  would  be  in  a 
similar  format  as  that  developed  for  the  airplanes.  Rearrangement  of  the  equation  results  in 
the  formulation  of  the  ground  friction  coefficient.  Vehicle  use  will  require  the  measurement 
of  stopping  distance  and  brake  application  speed.  Combination  of  this  data  with  all  other 
required  parameters  will  result  in  the  calculation  of  the  dynamic  ground  friction  coefficient. 
The  criteria  are  established  in  reference  to  the  developed  airplane  prediction  equation  and 
the  dynamic  process  it  describes.  Other  desirable  characteristics  for  such  a ground  vehicle 
are  ease  of  maintenance,  safe  and  simple  operation,  and  a data  collection  system  with 
minimum  complexity.  The  last  of  these  requirements  appears  to  be  met  by  existing  ground 
vehicles  such  as  the  Diagonally  Braked  Vehicle  (l)BV)  and  the  Mu-Meter.  However,  dynamic 
similarity  of  the  tires  used  is  absent.  The  DBV  is  a standard  automobile  that  allows  the  full 
lockup  of  one  front  and  one  rear  wheel.  Stopping  distances  are  measured  and  compared  for 
different  runway  conditions.  The  Mu-Meter  is  a trailer  that  derives  the  friction  coefficient 
through  a side  force  measurement  between  two  yawed  tires.  The  vehicle  is  towed  at 
constant,  speed  and  the  measurement  of  friction  coefficient  is  continuous.  Neither  vehicle 
duplicates  the  requirements  implied  by  the  prediction  equation.  More  details  and  a dimen- 
sional analysis  of  both  ground  vehicles  is  given  in  Appendix  F. 
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SECTION  XV 


EVALUATION  OF  DBV  AND  MU-METER 


For  the  prediction  equations  to  be  operationally  meaningful,  it  i necessary  to  accurately 
measure  ground  friction  (available  mu).  Two  ground  vehicles-NASA  Diagonally  Braked 
Vehicle  (DBV)  and  British  Mu-Meter  have  been  tested  extensive) . with  varying  degrees  of 
aircraft-to-vehicle  correlation  accuracy  (Refs.  2 through  8).  Both  vehicles  have  qualitative 
merit,  but  neither  has  the  quantitative  correlation  capability  r<  tuired  for  the  prediction 
equations. 

To  close  t.ie  gap  between  the  results  of  ground  vehicles  and  aircraft  measurements,  the 
physical  and  operational  reasons  for  the  differences  were  investigated.  Existing  ground 
vehicles  use  relatively  low  pressure,  light  motor  vehicle  tires,  while  all  aircraft  tires  are  high 
pressure.  Figure  28  and  Table  17  compare  an  aircraft  tire  and  a light  motor  vehicle  tire.  To 
expect  a good  correlation  between  the  two  tire-vehicle  combinations  without  accounting  for 
all  the  differences  shown  is  not  prudent.  Thus,  no  existing  ground  vehicle  meets  the  most 
fundamental  and  important  requirement  established  herein  of  using  a proper  tire. 


By  applying  model  laws  to  the  pi  terms  developed  earlier,  cert-i 
can  be  identified.  Subscripts  m denote  quantities  relating  to 
denote  quantities  relating  to  the  prototype.  Sufficient  condition 
Up  are  (see  Appendix  F): 

model  design  conditions 
odel,  and  subscripts  p 
he  equality  of  pm  and 

<cL,cD>m  = (tL/CD)p 

(59) 

vm  = <X',/6vp 

(60) 

sm  = (X)6  sp 

(61) 

where:  > = Wm/Wp,  which  is  the  weight  scale  ratio. 

For  reasons  explained  in  Appendix  F,  no  ground  vehicle  can  satisfy  Eq  59.  When  such  is  the 
case,  that  is,  when  a model  fails  to  satisfy  one  or  more  of  the  design  conditions,  it  becomes  a 
distorted  model.  This  can  be  properly  accounted  for  by  introducing  the  so-called  prediction 
and  distortion  factors  into  the  prediction  equation.  For  example,  for  the  current  problem, 
the  distortion  factor  is  of  the  form  (see  Appendix  F for  detailed  analysis): 

. -1C,-C4%SP. 

6 = a J * 

(62) 
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Table  17. -Comparison  of  Aircraft  and  Automotive  Tires 


Parameter 

Airplane  tire 

Ground  vehicle  tire 

Remarks 

Tire  foot- 
print 

Elliptical 

Rectangular 

Aircraft  tire  is  primarily  a torodial 
carcass  with  small  additional  tread 
rubber.  Automotive  tires  use 
heavy  tread,  particularly  in  the 
shoulder  region  (Ref.  9). 

Tire  thread 
design 

Natural  rubber 
base 

Synthetic  rubber  base 

Synthetic  rubber  base  materials 
can  produce  significantly  higher 
H on  wet  runways  (Ref.  10). 

Tire  tread 
design 

• Light  tread 

• High  stiffness 

• High  aspect  ratio 

• Heavy  tread 

• Low  stiffness 

• Low  aspect  ratio 

High  mileage  life  and  absorption 
of  road  surface  shocks  are  criteria 
for  automotive  tire.  Prevention  of 
excessive  loading  of  the  runway 
surface  is  major  criterion  for  air- 
craft tires  (Ref.  11). 

Tire  deflec- 
tion range 

30  to  35% 

12% 

Deflection  is  defined  as  percentage 
of  the  height  of  the  crown  above 
the  wheel  rim  (Ref.  11). 

Hydroplaning 
inception 
speed  (Vcr) 

•vcr=9y^r 

• Reversion  pos- 
sible for  all  tires 

•Vcr*  Kl/"psi 

for  constant 
deflections, 

• Reversion  possible 
only  if  pi  >50  psi  and 
v >60  mph 

Tire  stiffness  can  be  a factor  in 
determing  the  hydroplaning  in- 
ception speeds  for  smeller  size 
(model)  tires  (Ref.  12). 
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where: 


6 = prediction  factor 
a = distortion  factor 

and  C3,  C4  are  exponents  determined  experimentally.  See  Table  F-l  for  a list  of  distortion 
factors  to  be  used  with  various  prediction  equations  derived. 

When  these  criteria  are  applied  to  the  NASA  DBV,  it  fails  to  satisfy  Eq  5b,  and  no 
compensation  was  made  during  correlation  analysis  for  the  model  distortion.  When  the 
criteria  in  Eqs  60  and  61  are  applied  to  the  baseline  conditions  of  Tasks  I and  11,  the  results 
are  as  shown  in  Table  1 8. 

Table  18  shows  that  with  vehicle  weight  held  constant,  a different  velocity  is  necessary  each 
time  a different  condition  is  being  run.  Similarly  vm  could  be  held  constant  and  Wm  would 
become  a variable.  In  the  case  of  DBV,  however,  both  the  weight  and  velocity  were 
constants  at  5500  lb  and  60  mph,  thus  violating  model  laws.  This  situation  could  be 
corrected  by  running  the  vehicle  tests  for  a range  of  speeds  at  constant  weights  and 
establishing  a relationship  between  vehicle  speed  and  distance  under  various  runway 
conditions. 

A slightly  different  approach  was  tried  in  evaluating  the  mu-meter  to  apply  the  model 
theory  to  its  data.  Only  the  NASA  Wallops  airport  data  of  Ref.  3 was  evaluated.  Table  F-2 
shows  the  raw  data  taken  from  Ref.  3.  Since  the  data  shows  the  mu-meter  average  friction 
reading  for  each  run  and  the  corresponding  stopping  distance  used  by  the  727,  it  was 
decided  to  calculate  the  average  friction  (g)  used  by  the  airplane  by  applying  the  prediction 
equation  for  the  model  727  arrived  at  earlier  so  that  the  two  friction  readings  could  be 
compared.  The  details  of  calculating  ir  terms  are  shown  in  Table  F-3,  and  the  comparison  is 
shown  in  Table  19  and  plotted  in  Figure  29. 

The  comparison  clearly  indicates  the  randomness  of  the  data  and  the  missing 
quantitativeness.  Another  interesting  point  is  observed  by  comparing  runs  13  and  17.  Both 
runs  were  conducted  on  section  A-G  of  the  runway,  run  13  under  wet  conditions  (water 
0.01  ir..  deep)  and  run  17  under  dry  conditions.  The  mu-meter  indicated  a higher  friction 
reading  for  run  13  than  for  17,  while  the  resulting  airplane  tire  mu  was  closer  to  anticipated 
values. 

The  mu-meter  is  a continuous-recording  trailer  that  measures  the  side  force  generated 
between  the  test  surface  and  two  pneumatic  tires  that  are  set  at  a fixed  toe-out  angle  of  7.5 
degrees  to  the  line  of  drag.  The  available  side  load  thus  recorded  is  reported  as  a measure  of 
the  surface  friction. 

The  above  statement  implies  a few  broad  assumptions  that  may  not  be  true  for  actual 
aircraft  tire  operation.  For  example: 

• At  a 7.5-degree  yaw  angle,  the  tire  may  not  have  reached  the  maximum  and  limiting 
tire-runway  friction  condition  on  a dry  runway. 
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Table  18.-DBV  Correlation  Requirements 


Airplane 

Vp  (mph) 

Sp  (ft) 

Vm  (mph)' 

Sm  (ft) 

727 

133.0 

1278 

78.3 

443 

737 

118.0 

1068 

74.0 

421 

747 

149.3 

1805 

69.5 

413 

C-141 

136.4 

1841 

71.0 

500 

F-4 

174.6 

2768 

127.0 

1465 

*Wm  - 6200  lb  - WDBV 


19 


Aircraft 


1.0 


• Experimental  data  on  tire  cornering  on  wet  runways  show  that  the  cornering 
coefficient  reaches  a maximum  value  at  a critical  yaw  angle  and  then  decreases  with 
further  increase  in  yaw  angle.  Thus,  if  the  critical  yaw  angle  is  different  than  7.5 
degrees,  the  mu-meter  reading  will  underestimate  peak  a cornering. 

Ill  summary,  Figures  28  and  29  and  Tables  17  and  18,  and  19  clearly  show  the  inadequacy 
of  the  existing  ground  vehicles  in  meeting  the  previously  developed  criteria. 
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SECTION  XVI 

CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

The  following  conclusions  have  been  arrived  at  from  the  analysis  of  the  available  data 

gathered  during  this  program. 

• A dimensional  analysis  technique  can  successfully  express  braking  phenomena  in  the 
form  of  a mathematical  (model)  equation. 

• Experimental  data  from  an  airplane  braking  distance  sensitivity  study  are  needed  to 
determine  the  constants  and  exponents  in  the  model  equation  leading  to  a prediction 
model  equation. 

• With  the  proper  information  of  airplane  and  weather  parameters  and  an  accurate  and 
meaningful  measurement  of  the  tire-runway  friction  coefficient,  airplane  braking 
distances  can  be  predicted  within  reasonable  tolerances. 

• The  most  important  requirements  for  a vehicle  to  accurately  measure  tire-runway 
interface  friction  are  the  selection  of  a proper  tire,  (preferably  an  aircraft  type  of  tire) 
and  a faithful  reproduction  of  the  interface  dynamics,  such  as  the  use  of  a skid  control 
system. 

• Existing  ground  vehicles  (DBV  and  mu-meter)  fail  to  meet  the  criteria  described  herein 
and  are  distorted  or  dissimilar  models. 


2.  RECOMMENDATIONS 

If  a suitable  ground  vehicle  were  available,  a method  to  predict  stopping  distance  could  be 
outlined  as  shown  in  Figure  30.  The  ground  vehicle  is  operated  on  the  runway  to  measure 
brake  application  speed  and  stopping  distance.  Based  on  an  earlier  prediction  equation,  a 
ground  friction  coefficient  can  be  calculated  using  the  vehicle  parameters  and  runway  test 
data  as  an  input.  It  is  anticipated  that  the  resulting  ground  friction  coefficient  cannot  be 
directly  applied  to  the  airplane  prediction  equation  unless  some  data  on  tire  correlation  is 
applied.  Because  the  size  and  power  of  a ground  vehicle  are  limited,  the  test  tire  should  be  a 
small  aircraft  tire.  Because  aircraft  tires  of  many  sizes  are  being  used,  a scaling  factor  must 
be  applied  to  relate  aircraft  tire  performance  from  the  size  on  the  test  vehicle  to  other  sizes. 
It  appears  that  this  body  of  data  must  stili  be  generated. 

Weather  conditions  can  change  rapidly,  and  a repeated  and  frequent  prediction  of  the 
tire-runway  friction  coefficient  may  not  always  be  practical.  For  instance,  a short  but  heavy 
thundershower  may  change  a dry  runway  to  a flooded  one.  A correction  factor  for  water 
accumulation  on  the  runway  would  eliminate  the  need  for  an  immediate  new  ground  vehicle 
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Figure  30.- Recommended  Prediction  Methodology 
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friction  measurement.  It  is  anticipated  that  in  some  circumstances,  such  as  snow,  ice  and 
slush,  a weather  correction  factor  could  be  obtained  only  through  a direct  ground  friction 
measurement  by  the  vehicle.  The  runway  friction  coefficient  could  then  be  used  in  the 
prediction  equation  along  with  the  specific  airplane  parameters  to  calculate  the  stopping 
distance.  By  a comparison  to  available  runway  lengths,  a safety  margin  and  brake 
application  point  on  the  runway  can  be  established.  Finally,  the  prediction  model  approach 
needs  to  be  verified  on  additional  aircraft,  especially  fighter  types,  to  extend  the  application 
base. 
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APPENDIX  A 


CALCULATION  OF  APPROACH,  FLARE,  AND 
TRANSITION  DISTANCES 


Approach,  flare  and  transition  distances  were  calculated  with  the  aid  of  a digital  program. 
The  equations  used  in  the  program  along  with  the  assumptions  required  are  listed  in  suc- 
ceeding paragraphs.  Figure  A-l  helps  to  define  the  significant  variables  involved. 

1 . Approach  Distance 

Approach  distance  was  calculated  as: 


SA  = 57.30H/7 


where: 

SA  = approach  distance  in  feet 
H = height  above  threshold  in  feet 
7 = glide  slope  in  degrees. 

This  equation  assumes  no  speed  loss  during  the  approach  and  that  7 is  a small  angle  (for 
example,  tan  7 = 7). 

2.  Flare  Distance 

To  calculate  flare  distance,  a constant  centripetal  acceleration  has  been  assumed  to  exist 
during  the  flare.  This  assumption  results  in: 

SF  = 0.000271  V27/n-l 

where: 

SF 
V 
7 
n 


= flare  distance  in  feet 
= approach  speed  in  feet 
= glide  slope  in  degrees 
= flare  load  factor  > 1 .0 


per  second 

_ centripetal  acceleration 
gravitational  acceleration 
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3.  Transition  Distance 


During  the  transition  segment,  the  pilot  may  initiate  various  actions  that  can  affect  the 
aerodynamics  of  the  aircraft.  To  account  for  such  actions,  the  digital  program  used  a 
time-dependent  step  integration  routine  to  calculate  transition  distance.  The  basic  procedure 
involves  two  steps.  During  any  one  time  period,  an  aircraft  deceleration  is  first  calculated. 
Based  on  this  deceleration  and  the  initial  aircraft  speed  at  the  start  of  the  time  period,  a 
final  aircraft  speed  is  calculated.  The  distance  i«  then  computed  from  the  initial  and  final 
velocities  over  the  known  time  interval.  This  procedure  is  repeated  until  braking  is  initiated. 
The  applicable  equations  are: 

a = (FF.  - FD  - FB)/M 
VF  = VI  -(a  ■ T) 

DS  = 0.5  (Vi  +VF)T 
ST  *»  1 DS 

where: 

a = aircraft  deceleration 

FE  « aircraft  thrust 

FD  = aircraft  drag  force 

FB  = aircraft  braking  (rolling  friction)  force 

M = aircraft  mass 

T = time  interval 

VI  - initial  aircraft  velocity 

VF  = final  aircraft  velocity 

DS  “ distance  traveled  during  a time  interval 

ST  = total  transition  distance 

Airplane  thrust,  drag,  and  braking  forces  are  a function  of  vehicle  velocity.  Thus,  decelera- 
tion changes  throughout  the  transition  segment. 


APPENDIX  B 

BOEING  BRAKE  CONTROL  SIMULATION 


The  simulator  is  an  unalog-hardware  system.  The  aircraft  dynamics  arc  simulated  using 
analog  computer  equipment.  Wherever  possible,  actual  aircraft  components  are  used  instead 
of  simulation.  This  approach  is  used  to  provide  a more  accurate  simulator,  particularly  in 
areas  where  nonlinearities  and  complicated  dynamics  exist.  An  added  benefit  of  reduced 
computer  requirements  is  also  realized, 

The  simulation  can  be  divided  into  six  sections: 

• Hardware 

• Airplane  dynamics 

• Strut  dynamics 

• Wheel  dynamics 

• Brake  torque  simulation 

• Ground  force  simulation 

The  relationships  and  interactions  between  these  elements  are  shown  in  a detailed  block 
diagram  of  the  simulation  (Figure  B-l ).  The  detailed  mathematical  models  are  developed  in 
Ref.  1. 

1.  Hardware 

The  hardware  portion  of  the  simulation  consists  of  the  brake  system  hydraulics,  skid  control 
card  or  box,  skid  control  valve,  instrumentation,  and  signal  conditioning  equipment.  The 
hydraulic  portion  of  the  simulator  includes  all  hydraulic  components  that  can  influence 
brake  system  performance.  The  actual  hydraulic  system  mockup  typically  includes  the 
pilot's  metering  valve,  antiskid  valves,  accumulators,  pumps,  supply  and  return  lines,  reser- 
voir, and  the  brakes.  In  this  way,  the  hydraulic  system  characteristics,  including  any  non- 
linearity, are  accurately  simulated.  Actual  antiskid  control  circuit  cards  are  used  in  the 
simulator  to  ensure  that  the  overall  function,  nonlinearities,  component  characteristics,  and 
system  tolerances  are  accurately  reproduced.  Unless  these  components  are  reproduced  pre- 
cisely, the  usefulness  of  the  results  can  be  seriously  jeopardized.  The  signal  conditioning 
portion  requires  that  a DC  voltage  proportional  to  wheel  speed  (from  the  computer)  be 
converted  to  a frequency-modulated  signal.  The  FM  signal  is  the  input  to  the  control  box. 
Brake  pressure  is  monitored  and  becomes  an  input  to  the  computer  for  use  in  the  brake 
torque  simulation. 
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2.  Airplane  Dynamics 

A three-degree-of-freedom  model  is  used  to  simulate  airplane  dynamics.  The  model  includes 
the  vertical  and  pitching  degrees  of  freedom  in  addition  to  the  horizontal  degree  of  freedom 
(braked  airplane).  The  vertical  and  pitching  modes  are  used  to  account  for  dynamic  changes 
in  the  vertical  reaction  between  the  tire  and  runway.  The  braked  airplane  simulation  in- 
cludes the  effects  of  aerodynamic  forces,  engine  thrust,  and  ground  force. 

3.  Strut  Dynamics 

The  strut,  a very  important  aspect  of  the  simulation,  is  modeled  with  a two  degree  of 
freedom  system.  The  model  accounts  for  the  fore  and  aft  deflection  of  the  axle  centerline 
and  the  torsional  deflection  of  the  end  of  the  strut.  The  effects  of  the  strut  motion  on  the 
ground  force,  airplane  dynamics  and  measured  wheel  velocity  are  also  included. 

Special  attention  is  given  to  modeling  the  natural  frequency  and  the  damping  of  the  gear. 
This  is  done  by  determining  an  effective  strut  mass  and  fore-aft  stiffness  which  yields  a 
natural  gear  frequency  which  is  compatible  with  airplane  records. 

4.  Wheel  Dynamics 

The  wheel  dynamics  simulation  uses  the  ground  force  and  the  brake  torque  to  develop  a 
wheel  velocity.  The  simulation  of  ground  force  and  brake  torque  are  carried  out  by  separate 
models  as  described  in  the  paragraphs  below.  Data  used  in  this  portion  of  the  simulation  of 
the  tire  and  wheel  was  obtained  directly  from  the  tire  manufacturers'  published  data. 

5.  Brake  Torque  Simulation 

Although  the  actual  brake  is  used  for  the  hydraulic  load,  the  brake  does  not  serve  as  a 
torque  producing  element  in  the  simulator.  The  signal  generated  by  the  brake  pressure 
transducer  is  used  as  an  input  signal  to  the  brake  torque  simulation  which  includes  torque 
gain,  dynamic  response  retractor  spring  deadband,  torque  peaking  and  brake  fade.  The 
model  used  in  the  brake  simulation  is  the  result  of  extensive  research  including  both  full 
scale  brake  testing  (dynamometer  and  airplane)  and  small  scale  brake  lining  evaluation. 
Some  of  the  dynamometer  brake  tests  were  carried  out  under  two  recent  funded  Air  Force 
contracts  (F336 1 5-73-C-30 17  and  F3361 5-72-C-IOI5). 

6.  Ground  Force  Simulation 

The  ground  force  simulation  is  used  to  develop  the  horizontal  force  acting  between  the  tire 
and  ground.  The  method  used  is  based  on  an  experimentally  determine  ' relationship  be- 
tween the  ground  friction  coefficient  and  the  percent  slip  of  the  tire,  o.  Tire  heating  effects 
are  also  considered  in  the  model. 
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The  effects  of  a wet  or  contaminated  runway  are  also  reproduced  in  the  wheel  and  tire / 
ground  force  portion  of  the  simulation.  The  significance  of  vehicle  velocity  on  the  friction 
coefficient  and  the  path  clearing  effect  of  the  forward  tire  of  truck  type  landing  gear  are 
considered.  Data  from  Boeing  flight  tests  and  from  NASA/AF/FAA  tests  was  used  in  devel- 
oping this  portion  of  the  simulation. 

7.  Nonlinear  Relationships 

Several  nonlinear  relationships  were  used  in  the  aircraft  simulation.  These  relationships 
required  special  implementation  on  the  computer.  The  nonlinear  functions  used  were: 

• Wet  runway  curve 

• Touchdown  profile 

• Brake  torque-pressure  relation  for  two  conditions: 

• T = f(p)0'5 

• T = f(p)°'7J 

• Mu-slip  curves 

• Power  spectral  density  curve  of  rough  runway  test 

The  actual  relationships  used  are  documented  in  Figures  B-2  througu  B-6.  In  addition,  each 
relationship  is  briefly  described  in  the  following  paragraphs. 

la)  Wet  Runway  Curve.  - The  friction  coefficient  obtained  between  the  tire  and  ground  on  a 
wet  runway  is  a function  of  aircraft  velocity  (Figure  B-2).  Characteristically,  friction  in- 
creases as  airplane  speed  decreases.  The  wet  coefficient  of  friction  was  produced  on  the 
simulator  by  subtracting  a velocity-dependent  function  (Figure  B-2)  from  the  maximum 
value  of  friction  available  on  the  wet  surface. 

lb)  Touchdown  Profile. -The  touchdown  profile  as  implemented  on  the  simulator  is  a 
time-dependent  function  that  modifies  the  main  gear  load  (Figure  B-3).  The  function  was 
modeled  after  typical  bounce  conditions  seen  during  commercial  aircraft  landings. 

(c)  Nonlinear  Torque  Gain  Curves.  - The  static  torque-pressure  curve  of  a typical  aircraft 
brake  is  a nonlinear  function  (Figure  B-4).  The  curve  shape  is  dependent  on  brake  design 
and  the  physical  properties  of  the  brake  lining.  The  actual  curves  used  were  based  on  data 
produced  during  dynamometer  brake  tests. 

(d)  Mu-Slip  Curve  Variations.- The  baseline  mu-slip  curve  used  during  the  sensitivity  study 
represents  the  curve  commonly  accepted  throughout  the  aircraft  industry  (Figure  B-S).  The 
four  other  curves  were  derived  to  test  antiskid  system  reaction  to  changes  in  the  mu-slip 
relationship. 
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(e)  Power  Spectral  Density  Curve  of  Rough  Runway.- The  power  spectral  density  (PSD) 
curve  of  the  rough  runway  test  condition  is  a statistical  description  of  the  Kennedy  Airport 
runway  surface  (Figure  B-6).  The  area  beneath  the  PSD  curve  is  the  square  of  the  rms  value 
of  runway  roughness,  where  roughness  is  a measure  of  the  runway  vertical  profile.  The  shape 
of  the  PSD  curve  indicates  the  wavelength  of  a particular  roughness. 

8.  Recorded  Data 

In  addition  to  the  performance  indices,  the  following  data  was  recorded  on  an  eight-channel 
pen  recorder: 

• Braked  wheel  speed  (two  wheels) 

• Brake  pressure  (two  brakes) 

• Strut  displacement 

• Valve  signal  (two  wheels) 

The  variables  are  required  to  properly  assess  system  capabilities.  For  various  teats,  the 
following  weie  recorded  to  facilitate  analysis: 

• Supply  pressure 

• Return  pressure 

• Metered  pressure 

• Brake  torque 

• Ground  force 

• Vertical  displacement  of  airplane  CG 

The  ground  force  model  was  monitored  on  the  oscilloscope  to  observe  the  control  system 
operation.  This  was  necessary  to  initiate  some  stability  tests. 

9.  Test  Equipment 

A variety  of  test  equipment  and  recording  devices  were  required  during  the  study.  The 
equipment  used  throughout  the  testing  included: 

• Brush  chart  recorder 

• Digital  voltmeter 

• Weston  Boon  shaft  function  analyui 
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• Oscilloscope 

• X-Y  plotter 

• Wavetek  function  generator 

• Pressure  transducers  and  associated  devices 

Each  piece  of  test  equipment  was  calibrated  before  testing  began  and  as  required  thereafter. 
The  brush  chart  recorder  was  used  for  indication  only  and  was  not  calibrated. 

The  pressure  transducers  are  identified  in  Table  B-l. 


Brake  torque  {right  wheel) 


Figure  B-h- Block  Diagram  of  Boeing  Brake 
Control  Simulator 
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Of  Peak  Available  Mu 


Table  B-1  .—Pressure  Transducer  Information 


APPENDIX  C 

DIMENSIONAL  ANALYSIS  TECHNIQUE 


Dimensional  analysis  differs  from  other  types  of  analysis  in  that  it  is  based  solely  on  the 
relationships  that  must  exist  among  the  pertinent  variables  because  of  their  dimensions 
rather  than  the  laws  of  physics  and  classical  mathematical  derivation.  In  itself,  dimensional 
analysis  gives  qualitative  rather  than  quantitative  relationships,  but  when  combined  with 
experimental  procedures,  it  may  be  made  to  supply  quantitative  results  and  accurate  pre- 
diction equations. 

I . Form  of  Dimensional  Equations 

In  general,  any  measurable  quantity  a (a  secondary  quantity)  may  be  expressed  in  terms  of 
the  appropriate  quantities  a.  (called  primary  quantities)  that  affect  the  magnitude  of  the 
secondary  quantity.  The  general  relationship  between  a and  the  primary  quantities  may  be 
written  as: 


a = f fa | , an,  aj, an)  (C-l) 

in  which  a is  the  number  denoting  the  magnitude  of  the  secondary  quantity.  Therefore,  a| , 
a2,  a3 an  are  the  numbers  denoting  the  magnitude  of  the  significant  primary  quanti- 

ties involved. 

The  general  solution  (nature  of  the  function!  from  Eq  C-)  can  be  shown  to  be: 

“ = Caa|C‘  a2°2  a3C3 “nCn  (C'2> 

in  which  the  secondary  quantity  is  expressed  as  a dimensionless  coefficient  (Ca)  multiplied 
by  the  product  of  the  pertinent  primary  quantities,  each  raised  to  the  appropriate  power. 
Thus,  any  measurable  phenomenon  can  be  evaluated  by  Eq  C-2  in  terms  of  the  factors 
causing  it.  The  nature  of  coefficient  Ca  must  be  determined  experimentally. 

2.  Development  of  Prediction  Equations 

Development  of  prediction  equations  and  related  material  has  been  strictly  based  on  proce- 
dures described  in  Ref.  13.  As  a result,  Ref.  1 3 material  has  been  freely  used  and  quoted  as 
such. 

Two  general  methods  are  available  for  developing  prediction  equations.  One  method  consists 
of  establishing,  by  careful  observation  and  measurement,  the  effect  of  the  pertinent  varia- 
bles upon  the  quantity  to  be  predicted.  The  other  method  consists  of  applying  the  natural 
laws  pertinent  to  the  problem  to  develop  relationships  among  the  significant  variables.  The 
natural  laws  used  in  this  method  are  simply  generalizations  of  reliable  information  assem- 
bled through  observation  and  measurement.  The  first  method  is  usually  referred  to  as  the 
experimental  method,  the  second  as  the  analytical  method.  However,  each  involves  analysis, 
and  each  is  basically  dependent  upon  experimental  findings.  Often,  it  is  necessary  to  develop 
prediction  equations  for  phenomena  to  which  the  usual  analytical  procedures  are  not  well 
adapted. 
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The  first  step  of  determining  the  pertinent  variables  is  described  in  Section  X.  The  second 
step  is  to  express  the  secondary  quantity  as  a function  of  the  primary  quantities,  that  is: 

s = f (g,  v,  p,  Fe,  p,CL/CD)  (a) 

Next,  a relationship  is  developed  among  the  variables  by  dimensional  analysis,  leading  to  an 
expression  of  the  form: 

s « CgC>  v*"2  F^4  pCs  (CL/CD1  6 (b) 

The  determination  of  the  prediction  equation  for  s,  from  Eq  b,  would  involve  the  evaluation 
of  seven  unknowns  (C,  C(,  C2,  C3,  C4,  C5,  and  Cg)  but  the  number  of  unknowns  can  be 
reduced  with  the  aid  of  dimensional  analysis.  The  dimensional  equation  corresponding  to  Eq 

The  variables  p and  C l/C q are  dimensionless.  Equation  c can  be  resolved  into  three  com- 
ponent auxiliary  equations: 


for  M . 

0 s*  C3  + C4 

(d) 

for  L: 

1 = C,+C2-3C3  + C4 

(e) 

forT: 

0 « -2C,  -c2-  2C4 

O') 

There  are  now  four  unknowns  and  three  equations:  hence  three  of  the  unknowns  may  be 
expressed  in  terms  of  the  one  remaining  unknown.  Several  combinations  are  possible;  one 
such  combination  is  to  express  C 1 , C2,  and  C4  in  terms  of  C3.  From  Eq  d,  e,  and  f: 


C,  = -1-2C3 
C-,  = 2 + 6C3 
c4  = -C3 

It’  these*  values  are  substituted  into  Eq  b,  the  result  is: 

,.C  (A  }(,%<■>) 

This  may  be  consolidated  to: 


(g> 

(It) 

(i) 


(j) 


(k) 
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The  original  equation  involving  seven  unknowns  has  been  reduced  to  an  equation  involving 
only  four  unknowns,  which  represents  considerable  simplification.  Several  forms  of  the 
reduced  equation  are  possible,  dependent  on  which  exponents  are  retained.  Other  possible 
forms  are: 


.-qk5l  m 


(CL/cd)C 


MCs  (cl/cd)C 


In  each  case,  the  quantities  in  parentheses  are  dimensionless.  The  coefficients  Ca,  Cp,  and 
Cy  are  dimensionless  and  are  functions  of  the  dimensionless  groups: 


ca=  (*)  /JL«,  ,p,cL/cD 


From  this  it  follows  that  Eq  k can  be  written: 


(5)-  ■■  f • 


Equation  C-4  is  a simplified  form,  involving  an  unknown  function  of  three  variables  instead 
of  the  seven  unknowns  of  Eq  b.  Each  term  in  Eq  C-4  is  dimensionless. 


Thus  the  term  CQ  is  seen  to  be  a function  of  dimensionless  groups  of  the  variables  influ- 
encing the  phenomenon.  The  dimensionless  groups,  which  are  known  as  pi  terms,  may  be 
designated  as  iri.  Hence,  Eq  C-4  can  be  written  in  the  general  terms: 


ir,  = F (tr2,ir3,  • 


in  which  k denotes  the  total  number  of  dimensionless  groups  affecting  the  phenomenon. 
The  number  of  dimensionless  groups,  or  pi  terms,  required  to  express  a phenomenon  can  be 
determined  from  the  Buckingham  pi  theorem. 


J 


;T  ?'7T|WA>W' 


■z  -( •.:,^VT>'- ' T-'p'-'Z'rT.l  T'vpr-W’rsrvTVT r_*3  p p J „,T 
•••  « ' : 


3.  The  Buckingham  Pi  Theorem 

This  section  has  been  developed  based  on  Ref.  14  material.  The  Buckingham  pi  theorem, 
which  is  the  basis  of  nondimensional  analyses,  asserts  that  any  complete  physical  relation- 
ship can  be  expressed  in  terms  of  a set  of  independent  dimensionless  products  composed  of 
the  relevant  physical  parameters.  Bridgman's  more  mathematical  statement  of  the  pi  theo- 
rem says,  “If  the  equation  F (a] , 02,  03 an)  = 0 is  complete,  the  solution  has  the 

form  f (*|,  7T2 tfn-k)  = where  the  it  terms  are  independent  products  of  the 

parameters  a\ , 02,  etc.,  and  are  dimensionless  in  the  fundamental  dimensions.”  The  number 
of  pi  terms  in  the  solution  equation  is  less  than  the  number  of  parameters  by  a factor  k. 
Usually  k equals  the  number  of  fundamental  dimensions;  however,  under  certain  circum- 
stances, k can  be  less. 

4.  Statement  of  Dimensional  Homogeneity 

Because  pi  terms  are  products  or  quotients  of  the  original  parameters,  and  because  these 
products  or  quotients  are  of  zero  dimension,  we  can  write  a general  dimensional  equation 
expressing  this  fact.  The  most  general  grouping  of  the  parameters  entering  the  current 
problem,  for  example,  must  combine  in  such  a way  that  the  dimensions  of  the  product  are 
zero;  that  is: 


sCl  %2  vC3  pC4  Ft,Cs  b (Ci/x))C7  i M°  L°  T°  (C-6) 

d 

Equation  C-6  is  called  an  equation  of  dimensional  homogeneity.  The  symbol  = means 
"dimensionally  equal  to.”  The  symbols  M,  L,  and  T are  the  three  fundamental  dimensions  in 
the  physical  system  of  mass,  length,  and  time  respectively.  Equation  C-6  implies  that  the 
product  of  all  pi  terms,  each  of  which  is  singularly  nondimensional  of  its  own  accord,  will 
also  be  of  zero  dimension  when  expressed  in  the  fundamental  units  of  measure. 

S.  Determination  of  Pi  Terms 

The  easiest  and  most  orderly  procedure  for  developing  pi  terms  from  a list  of  variables  is  as 
follows:  The  seven  parameters  needed  to  define  the  problem  are  summarized  in  Table  C-l , 
where  their  fundamental  dimensions  are  also  presented  in  the  mass  (M),  length  (L),  and  time 
(T)  system.  The  list  of  parameters  to  be  manipulated  can  be  reduced  from  seven  to  five  by 
writing,  by  inspection,  two  pi  terms  that  are  already  nondimensional  quantities. 


7T|  = p and  it 2 = CL/CD 
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Next  we  arrange  the  remaining  five  quantities  in  a matrix  with  their  fundamental  dimen- 
sions. Across  the  top  cf  the  array  we  write  the  variables,  and  down  the  left  side  we  write  the 
fundamental  dimensions,  in  this  case  M,  L,  and  T.  Under  each  variable,  we  write  the  powers 
to  which  each  dimension  is  raised  in  each  variable. 


s g v p Fe 


M 0 0 0 11 

L 111-31 

T 0-2-10-2 


We  now  apply  the  matrix  algebra  theorem,  which  states  that  "from  a matrix  there  will  be  a 
number  of  independent  equations  equal  to  the  rank  of  the  matrix”;  the  rank  of  the  matrix  is 
defined  to  be  the  order  of  the  highest-order  determinant  of  the  matrix  that  differs  from 
zero.  Because  in  dimensional  analysis  there  always  exists  more  columns  than  rows  and 
columns,  we  can  select  numerous  combinations  of  different  columns  from  the  matrix  and 
combine  them  into  determinants.  Some  of  the  determinants  from  our  matrix  are: 


0 0 0 

0 0 1 

0 1 1 

0 1 

1 0 

1 -3 

1 1 1 

1 1 -3 

1 -3  1 

1 1 

-3  1 

0 0 

0 -2  -1 

-2  -1  0 

-1  0 -2 

(a)  (b)  (c) 


The  order  of  a determinant  is  the  number  of  columns  or  rows  in  the  determinant.  The  first 
three  determinants  are  third-order  determinants,  and  the  last  three  are  second-order  determi- 
nants. Although  the  determinant  of  (a)  equals  zero,  the  determinant  of  (b)  equals  + 1 .0  and 
is,  therefore,  not  zero.  This  observation  means  that  the  rank  of  our  matrix  is  3 and  that  the 
total  number  of  dimensionless  products  in  a complete  set  is  equal  to  the  total  number  of 
variables  minus  the  rank  of  their  dimensional  product  or  minus  the  number  of  independent 
equations.  Thus  we  have  5-3,  or  2 more  ir  terms.  Our  three  equations  can  be  written  by 
inspection  from  the  preceding  matrix  as: 


C4  + C5  = 0 

(C-7a) 

Cj  + C2  + C3  • 3Cq  + C5  = 0 

(C-7b) 

-2C2-C3-2C5  = 0 

(C-7  c) 
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Solving  in  terms  of  C2,  and  C3,  and  C 5 we  get: 


C2  = C,  - 2C4 
C3  = 6C4-2C, 

C5  = ‘C4 

Substitution  for  C2,  C3,  and  Cj  in  a statement  of  dimensional  homogeneity  gives: 


SC1  gCr2C4v6C4‘2Cl  pC4  Fe'C4  = MoLoTo 


Collection  of  exponents  with  the  same  coefficient  gives: 


C4  = M°L°T° 


These  two  nondimensional  quantities  are  the  additional  pi  terms.  Rearrangement  of  the  two 
pi  terms  obtained  algebraically,  together  with  the  two  pi  terms  written  by  inspection,  gives: 


IT|  = (sg/v2) 

Ir2  = M 
*3  = (cL/cd) 

77  4 = (pv6/Feg2) 

This  result  yields  the  functional  relationship  that: 


(sg/v2 ) = F (p,CL/CD,pv6/Feg2) 


This  equation  is  identical  in  form  with  Eq  C-4  solely  because  of  the  selection  of  the  group  of 
exponents  in  terms  of  which  the  independent  equations  were  solved.  Actually,  pi  terms  can 
be  multiplied  to  form  new  groups  of  pi  terms.  They  can  be  inverted ; they  can  be  squared,  or 
their  square  roots  can  be  taken.  Such  manipulations  are  proper  and  are  usually  performed  to 
create  a more  convenient  ratio  of  physical  phenomena.  Most  often,  one  isolates  the  de- 
pendent variable  in  a single  pi  term. 
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Table  C~1  .-Parameters  for  Braking  Stop  Distance  Model 


Fundamental  dimension 


Parameter 

Symbol 

Braking  distance 

s 

Acceleration  of  gravity 

g 

Brake  application  speed 

V 

Air  density 

p 

Forward  thrust 

% 

Coefficient  ot  friction 

p 

Coefficient  of  lift 
Coefficient  of  drag  ratl° 

cL/cD 

APPENDIX  D 

FORMULATION  OF  COMPONENT  EQUATIONS 


When  the  experimental  data  had  been  arranged  as  described  in  Appendix  C,  relationships 
between  irj,  the  term  containing  th<  dependent  variable,  and  ir 2,  713  and  77 4 in  turn,  the 
terms  with  independent  variables,  were  obtained  using  statistical  curve  fitting  programs.  The 
relationships  between  ir|  and  other  individual  it  terms  are  called  component  equations. 

Plots  w»re  prepared  of  tt|  vs  7T2,  tt  | vs  773,  and  tt  j vs  774  for  all  airplanes  using  data  from 
ASD-TR-74-41,  Volume  11,  Section  XIII.  An  example  of  these  plots  is  shown  (for  737  data) 
in  Figure  D-l.  This  helped  determine  the  general  form  of  relationship  that  could  exist 
between  vr  ] and  7T2,  7r  j and  7T3,  and  so  on;  e.g.  the  tt;  vs  773  data  plotted  as  a straight  line  on 
the  log-log  paper  and  therefore  should  have  a relationship  of  the  form 

y = AxB  (a) 

where  A is  a constant  and  B is  a polynomial.  Logarithms  of  both  sides  in  Eq  a give: 

!n(y)  = in  A + B In  x 


or: 


In  (77| ) = In  A + B In  (7^) 

which  is  the  equation  for  a straight  line,  Thus,  with  tt|  and  tt2  as  inputs  and  the  desired 
output  in  the  form  of  Eq  a,  two  computations  were  necessary,  namely  a log  transformation 
and  a determination  of  the  constants,  the  latter  requiring  polynomial  regression. 

1 . Polynomial  Regression 

Polynomial  regression  is  a statistical  technique  for  finding  the  coefficients  in  a functional 
relationship  between  a dependent  variable  (y)  and  a single  independent  variable  (x).  Powers 
of  an  independent  variable  are  generated  to  calculate  polynomials  of  successively  increasing 
degree.  If  there  is  no  reduction  in  the  residual  sum  of  squares  between  two  successive 
degrees  of  polynomials,  the  calculation  will  be  terminated  before  completing  the  analysis  for 
the  highest  degree  polynomial  specified.  This  essentially  is  a method  to  establish  minimum 
error  in  the  solution. 

Tlie  component  equations  for  TT]  vs  7^  and  7T]  vs  774  were  determined  in  the  fashion  just 
described.The  curvature  of  the  jt  1 vs  773  plots  (Figure  D-l)  indicates  that  a straight  line 
relationship  such  as  Eq  a would  not  yield  the  desired  accuracy.  Thus,  it  was  necessary  to  use 
multiple  linear  regression  in  addition  to  the  log  transformation. 


145 


2.  Multiple  Linear  Regression 


This  is  a statistical  technique  for  analyzing  a relationship  between  a dependent  variable  (y) 

and  a set  of  independent  variables  (xj,  x2 xn).  Least-squares  methods  are  used  to 

estimate  the  coefficients  in  a linear  relationship.  The  analysis  showed  that  the  relationship 
will  be  of  the  form: 


y = Ax,<b-Cx2> 

where  X2  was  found  to  be  % spoiler  (configuration). 

Thus  an  intermediate  calculation  was  necessary  between  log  transformation  and  multiple 
regression  to  determine  the  product  (spoiler  percentage)  In  (irj).  A flow  chart  (Figure  D-2) 
depicting  the  formulation  of  component  equations  follows  Figure  D-i . 

The  curve  fitting  computer  programs  were  used  from  BCS-STATPK,  Ref.  15. 
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igure  D- 1.  —Plots  of  it  y vs  tt^,  itj,  and  tt#  (737  Data) 


Figure  D -2. —Equation  Flow  Chart 


. | 


APPENDIX  E 

FORMULATION  OF  GENERALIZED  FUNCTIONS 

1 . Determination  of  Functions 

As  pointed  out  in  section  XI,  the  component  equations  representing  relationships  between 
ttj  and  the  other  individual  pi  terms  are  combined  to  give  a general  relationship.  This 
combination  is  not  always  simple,  but  under  certain  conditions  it  may  be  reasonably  direct. 
The  function  on  the  right-hand  side  of  the  general  equation: 


it  j - F (tt 2,  fl  ^4, 


■ V 


(E-l ) 


may  denote  any  combination  of  the  pi  terms.  We  will  now  determine  the  conditions  neces- 
sary or  sufficient  for  certain  simple  combinations  to  exist. 

2.  Conditions  for  Function  to  be  a Product 

If  four  pi  terms  are  involved  in  a phenomenon  (as  in  the  current  situation): 

tT|  = F (rr2,  irj, 7T4>  (E-2) 

experiments  would  be  carried  out  varying  Jr2  an(*  holding  and  constant.  From  a plot 
of  tr 1 against  rr2,  the  relationship: 


(IM)3,4  ~ fl  (,r2’  S3’ f4> 


(E-3) 


in  which  the  bar  denotes  constant  values,  could  be  established.  From  another  set  of  experi- 
ments, with  both  tr2  and  #4  constant  and  ^3  variable: 


(ir|)2,4  = f2  (tt2, 7-3, 7F4) 

may  be  established.  Similarly,  the  relationship: 

(ff  I ) 2,  3 “ *3  ^2'  *3’  *4) 


(E-3a) 


(E-3b) 


could  be  established.  Equations  such  as  E-3,  3-3a,  and  E-3b,  determined  by  holding  all  but 
one  of  the  pi  terms  in  the  function  constant,  are  called  component  equations. 

Under  certain  conditions,  the  component  equations  can  be  combined  to  form  the  general 
prediction  equation  by  multiplication.  For  example: 


*1  = C(’rl)3,4('rl)2,4(7rl*2,3 


(F.-3c) 


To  establish  those  conditions,  we  first  determine  the  constant  C in  Eq  E-3c  by  assuming  that 
the  component  equations  are  simply  multiplied  to  form  the  general  equation: 


! 4 


\ 1 


I i 
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F (tt2,  ff3,  W4)  = f(  (jt2.  W3,  ?4)  f2  <w2,  n-3>  ^4)  f3  (w2,  W3,  W4)  (E-4) 

If  this  is  true,  the  first  set  of  tests,  with  tr3  and  tr4  constant,  will  give: 

F (ir2, 5r3,  jF4)  = f,  (tr2,  ?3,  f4)  f2  (iF2,  ff3,  ?4)  f3  (?2,ff3,?4)  (E4a) 


from  which: 


F (tr2, 7r3,  tr4) 

fl  <'1’2'  ^4*  f2(ff2,tf3,  ff4)  f3(ir2,  ir3,  C;) 


The  second  set  of  tests,  with  nj  and  ir^  constant,  gives  from  Eq  E-4: 


F (7f2,  jr3,  if4>  = f)  (if2,  ff3,  tf4)f2  (f2,  tr3,  ff4)  f3  (tr2,  #3,  ?4)  (E4c) 


from  which: 


F (?2,  ir3,  f4) 

f2  (tt2,  w3,  v4)  (ir2,  f3,if4)  f3  (tf2,  *r3,  tr4) 


The  third  set  of  tests,  with  rr2  and  tr3  constant,  gives,  from  Eq  E4: 

F (ff2,ff3,tr4)  = f,  (f2, 5F3,  ff4)  f2  (tf2,#3,»4)  f3  (#2,f3,tr4)  (E4e) 


F (tr2,  sr3,  tt4) 

f3  (^2*^3>*4)  = 7 72  - - , c 7Z — z err 

J J ^ f|  (#2,  iT3,  w4)  f2  (tt2,  tt3,  tr4) 


values  of  f j (tr2,  j?3,  if4),  f2  (tt2,  tt3,  jf4),  and  f3  (ir2,  ?3,  ir4)  from  Eqs  e4b,  e4d  and  e4f 
are  substituted  into  Eq  E4  to  give: 


F (rr2,  jt3,  ?r4)  = 


F (tr2,  f3,  rr4)  F (jr2, 7r3,  f4)  F (#2,  jf3,  ir4) 
[fl  (ff2'  f3>  f4>  f2  (ff2'ff3’S4*  f3  (ff2'S3’ff4'] 


2 (E4g) 


?ti’:vrT*WFr^*?>K,v rs? y r~  '"wtyw?  - 


^2ii2SB!iSS££ii^a££2SSSiSS5SffiSI^ii 


8»>.  - ■-  *P;  "•'jaj 


However,  the  denominator  of  Eq  E-4g  is  found  from  Eq  E-4  with  all  773,  773,  and  7r4 
constant: 


F (jtt,  7F3,  ir4)  = f|  (f2.  S3, 14)  f2  (tf2,  S3,  S4)  f2  (tr2,  Sj.sq) 


(E-4h) 


Hence: 


F (773, 773,  774) 


_ F (7T2,  ^3,  ^4)  F (f2, 773,  S4)  F (ffo,  W3,  tt4) 
[F  (5f 2,  #3,  S4)]2 


(E-4i) 


In  addition  to  giving  the  value  of  C in  Eq  E-3e  as  l/|F(f2,  S3,  S4)]2  Eq  E-4i  indicates  that 
the  three  component  equations  must  have  the  same  form. 

A test  for  the  validity  of  combining  the  component  equations  as  a product  may  now  be 
developed  by  assuming  that  a fourth  component  equation  is  determined  from  a fourth  set  of 
data  in  which  one  of  the  pi  terms  is  held  constant  at  a different  value  than  in  the  preceding 
set  of  data.  For  example,  the  general  EqE-4i  was  determined  by  holding  tt2  constant  at  a 
value  of  ?•>,  but  if  valid,  it  could  also  have  been  determined  from  a set  of  data  in  which 
tr2  = a i . Then : 


„ F (tr-),  TTs,  tr4)  F ("•),  7M,774)  F (7?^,  Sq,  W4) 

F <ir2, 773,  tt4)  = - V - ' :3~2 

[F  (7T2,  7T3,  jt4)J  l 


(E-5) 


The  right-hand  side  of  Eq  E-4i  must  equal  the  right-hand  side  of  Eq  E-5.  Hence: 


F (S3, 773,  S4)  F (771,773,774!  F (773. 773. 774)  F (S2. 773, 774) 
[F  (if2.773. 774)] 2 [F  <*2, 773.  >74)]  - 

Similarly,  if  7r2  had  been  held  constant  at  a different  value,  773: 

F (77-1,773,774)  F (773,773,774)  F (773,773.774)  F (ffi.  773. 774) 
[F  (77-t,  7f  3,  Sql]  - [F  (S3.  773.  774)]  3 


(E-6) 


E-ha) 


Equations  E-6  and  E-ha  constitute  a test  for  the  validity  of  Eq  E-4i.  That  is,  if  the 
supplementary  sets  of  data  satisfy  either  Eq  E-6  or  E-ha,  the  general  equation  may  be  formed 
by  multiplying  the  component  equations  together  and  dividing  by  the  constant,  us  indicated 
in  Eq  E-4i. 

Thus,  if  the  general  equation  for  a system  involving  k pi  terms  formed  by  multiplication  of 
the  component  equations,  it  may  be  shown  that  the  form  is: 


77  1 = F (773,  773,  774 ff^)  F (773.773.7^3 77^)....  F ( 773,773,774..,,,  77^1 

[F  (f2, 773,  )]  k‘“ 


( E-7 ) 
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In  this,  tt2  must  be  the  same  in  each  component  set,  ^3  must  be  the  same  in  each  set,  etc. 
Similarly,  if  the  general  equation  for  a system  involving  k pi  terms  is  formed  by  addition  of 
the  component  equations,  it  may  be  shown  that  the  form  is: 


*1  = F (»2>53,?4 *k>  + F (»2,ff3’*4 *k)  + F (S2.Tr3,tr4 ti^)  {Fg) 

+ ....  F (ff2,ff3,)?4 nj( ) - (k-2)  F Uf2,n3,ir4 fy) 

In  general,  the  possible  methods  of  formation  of  genera!  equations  discussed  above  are 
adequate  for  the  majority  of  engineering  problems.  Regardless  of  whether  the  resultant 
prediction  equation  is  formed  by  multiplication  or  by  addition,  a constant  term  of  the 
form: 


Ffjf2.S3.if4 ”k> 

is  involved,  and  this  constant  can  be  evaluated  from  any  one  of  the  component  equations. 
As  a general  policy,  the  constant  should  be  evaluated  from  each  of  the  component  equa- 
tions. Each  should  give  the  same  value  as  the  others.  If  not,  error  is  present,  and  the 
equations  should  be  checked.  Tables  E-l  and  E-2  show  detailed  calculations  for  tests  of 
validity. 


Table  E-2.— Calculations  of  Validity  for  Constant  Term 


4.293  4.222  - 0.05647  0.05426  0.0561 


APPENDIX  F 

MODEL  THEORY  ANALYSIS  FOR  GROUND  VEHICLES 


1 . Model  Theory 

Two  coirmon  systematic  methods  can  derive  modeling  laws  for  a particular  physical  system. 
The  “equations  approach"  requires  knowledge  of  the  characteristic  equations  that  govern 
the  behavior  of  the  system.  The  “parameters  approach”  requires  only  that  the  complete  set 
of  all  variables  affecting  the  system  behavior  be  specified.  For  our  problem,  the  second 
approach  is  more  appropriate. 

We  have  already  considered  the  independent  variables  affecting  braking  stop  distances;  by 
applying  conventional  techniques  of  dimensional  analysis,  we  have  established  the  relation- 
ship: 

(sg/v2)  - F,  (it,  CL/CD,  Pv6/FeJ) 

However,  we  can  rewrite  the  above  equation  as: 

<**)  = f2  (cL/cD . sg/v2,  pv6/Fe  g2)  (p.,. 

The  model  design  conditions  follow  directly  from  Eq  F-l.  By  letting  the  subscript  m denote 
quantities  relating  to  the  model  and  the  subscript  p denote  quantities  relating  to  the  proto- 
type, we  find  that  sufficient  conditions  for  the  equality  of  and  are: 


Fem 


1 


The  conditions  given  by  Eqs  F-3  and  F-4  can  be  met  under  a variety  of  test  constraints. 
Equation  F-2,  however,  presents  a unique  problem.  No  ground  vehicle  can  duplicate  the 
wide  variety  of  Cj_  and  Cp  for  a single  airplane,  and  certainly  not  for  all  the  airplanes. 
Because  of  this  problem,  it  is  not  practically  feasible  to  conduct  a model  test  without  some 
measure  of  model  distortion.  A systematic  technique  to  account  for  model  distortion  has 
been  developed  by  Murphy  Ref.  13  and  will  be  used  later. 

The  model  can  be  tested  with  the  same  fluid  as  the  prototype  under  the  same  gravitational 
conditions.  That  is: 


(F-5) 

and 

Pm  = Pp 

(F-6) 

Simultaneous  solution 

of  Eqs  F-3  through  F-6  results  in 

the  following  set  of  simplified 

design  conditions: 

/ Fe  \1/6 

-v 

p 

(F-7) 

=(^t)sp 

(F-8) 

In  generating  Eq  F-l , if  we  also  considered  the  weight  of  the  vehicle  as  an  independent 
variable,  the  relationship  is  altered  as  follows: 


<M>  = F3 


fcL/CD.li : ,£L6 

\ v2  Wg2  W/ 


from  which: 


vm^  _ Ppvp6 

Wm  gm')  W„  B 2 
m p 


(F-4a) 


H 
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and: 


(F0  (\) 

"w "w — (F-4b) 

m "p 

Substituting  Eq  F-4b  into  Eq  F-4a  gives: 

Pm  ^ _ lpvp6 

Fem8m2"  Fepgp2  (F-40 

Equation  F-4c  is  identical  to  Eq  F-4.  Hence,  Eq  F-7  can  be  modified  as: 


v 


m 


v 


P 


or: 


vm  = Wl/6  vp  (F-9) 

and: 

sm  = <x>6  sp  (F-10) 

Where  X = Wm/ Wpis  the  weight  scale  ratio  tor  the  experiment. 

2.  Model  Distortion 

A distorted  model  is  one  in  which  one  or  more  of  the  design  conditions  are  not  satisfied. 
That  is,  one  of  the  pi  terms  is  not  equal  to  the  corresponding  pi  term  in  the  prototype.  In 
general,  the  prediction  equation  is  formulated  by  dividing  the  general  equation  for  the 
prototype  by  the  general  equation  for  ihe  model: 


*1  f (*  2’  * 3 . tr4) 

#lm  f (,r2nv  11 3w  Mm* 
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If  the  des  1 nditions  for  a true  model  are  all  satisfied  (ffjm  = irj),  the  functions  are  equal 
and  7T  j However,  if  one  design  condition  is  violated,  the  functions  may  not  be  equal, 

whereupon  it  j will  not  equal  7Tjm. 

Under  those  conditions,  two  procedures  are  available  for  establishing  a prediction  equation: 
• Determine  a prediction  factor  6 so  that: 


= 6trlm 

• Distort  additional  pi  terms  in  a controlled  fashion  so  that: 


(F-12) 


I 


'!m 


(F-13) 


(a)  The  Prediction  Factor,  8 - If  one  of  the  design  conditions  is  not  satisfied  (say 
irj),  the  prediction  factor  is,  by  definition  from  Eq  F-l  1 : 


^ "1  f (tr2,  ir3,  tr4) 

^lm  ^*2m*  *3m*  *4m) 


(F-l  2a) 


Hence,  to  evaluate  8,  the  ratio  of  the  two  functions  must  be  evaluated.  This  determination 
involves  either  additional  experimental  evidence  or  knowledge  of  how  irj  influences  the 
function.  If,  for  example,  it  can  be  established  that: 


f (»r,,  irv  ir4)  = 


f (tr^)  f (tr2,  tt4) 

f (TT,.  Tty  1 Tq) 


(F-l  4) 


it  follows  that: 


6 * 


r <*3> 

f <»3m> 


(F-l  2b) 


The  test  for  the  resolution  of  a function  into  the  product  of  two  component  functions 
requires  that  two  sets  of  tests  be  run  in  which  the  degree  of  distortion  is  varied.  The  degree 
of  distortion  can  be  evaluated  as  a distortion  factor  a,  defined  us: 


*3m  " an3 


(F-15) 


Thus  a = 1 indicates  no  distortion. 
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If  the  log  plot  of  7T]  vs  ir3  is  a straight  line  and  rr3  satisfies  the  requirement  for  combination 
as  a product,  it  follows  from  Eq.  F- 1 2b  that: 


f (1t3) 

6 - 

f (air3) 

(F-l  2c) 

C ir3m 

6 = 

(F-l  2d) 

Cam  „3m 

5 = a"m 

(F-l  2e) 

where  m is  the  slope  of  the  line  in  the  log  plot,  Even  if  values  of  7Tj  plotted  against  tr 3 do 
not  plot  as  a straight  line,  the  value  of  IT|,  corresponding  to  the  correct  value  of  W3  for  a 
true  model,  can  be  taken  from  the  curve,  if  a number  of  points  have  been  obtained  for 
values  of  jr3,  both  less  and  greater  than  the  desired  value,  and  if  the  curve  is  smooth,  the 
value  of  it  | may  be  determined  with  a high  degree  of  accuracy. 

(b)  Multiple  Distortion. -In  some  situations,  distortion  of  one  factor  or  dimension  will 
result  in  the  distortion  of  more  than  one  pi  term;  two  or  more  factors  may  be  distorted, 
causing  two  or  more  pi  terms  to  be  distorted.  The  prediction  factor  6 then  becomes  a 
function  of  the  distortion  factors,  and  may,  in  addition,  be  a function  of  one  or  more  of  the 
pi  terms.  It  may  be  determined  either  algebraically  or  experimentally. 

(cl  Compensated  Distortion  - One  simplification  of  the  prediction  equation  that  is  some- 
times possible  is  to  adjust  the  distortion  factors  in  such  a manner  that  the  prediction  factor 
becomes  unity.  Then: 


*1  ~ W1  m 


(F-13) 


as  it  would  be  in  an  undistorted  model. 

If  enough  is  known  concerning  the  effect  of  each  of  the  pi  terms,  this  can  be  accomplished 
algebraically  without  recourse  to  experiment  (see  Figure  F-l).  In  the  current  problem,  rr3  # 
ir3m.  Using  Eq  F-l  2e,  we  can  write: 


s =-|C3-C4%SP] 


(F-l  6) 


where: 


8 = f(ir3)/f(rr3m) 
a = Ir3nV’r3 
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Here,  C3  and  C4  are  the  constants  determined  in  the  preceding  prediction  equations.  Based 
on  Eq  F-16,  prediction  and  distortion  factors  corresponding  to  prediction  Eqs  4 through  38 
arc  listed  in  Table  F-l.  Table  F-2  shows  the  data  used  for  mu-meter  evaluation,  and  Table 
F-3  shows  its  detailed  analysis. 


Table  F-1.—  Prediction  Factor  3 vs  Distortion  Fictor  a 
for  DBV  Evaluation 


To  be  used 
with  Eq : 


5 » a( 0.337 5 % SP  * 0.3125] 

(4) 

j = a(0,301t  % SP  - 0.309] 

(3) 

5 . a(0.1347  %SP-  0.2159] 

(9) 

5 . a[0.081 13  % SP  -0.15668] 

(12) 

j = aI0, 08663  % SP  • 0.16666] 

(14) 

5 - a[0. 06279  % SP  • 0.1 6397] 

(16) 

j = a{0.1 1286  % SP  0.14803] 

(19) 

5 „ a[0.09193%SP  - 0.18967] 

(21) 

5 . a[0.081 8 %SP-  0.2365] 

123) 

5 = a(0. 13084  %SP  0.21341]  1 

126) 

5 = a[0.1 4819  %SP-  0.23065] 

(28) 

{ = a[0.16567  % SP  - 0.21614] 

130) 

6 = a'0-383® 

(33) 

6 = “[0.52698  + 0.08287  %SP) 

(36) 

6 = <,-0.89738 

Calculated 


